CARMICHAEL NUMBER VARIABLE RELATIONS: 
THREE-PRIME CARMICHAEL NUMBERS UP TO 10 



J.M.CHICK 

Abstract. Bounds and other relations involving variables connected with 
Carmichael numbers are reviewed and extended. Families of numbers or in- 
dividual numbers attaining or approaching these bounds are given. A new 
algorithm for finding three-prime Carmichael numbers is described, with its 
implementation up to 10^^. Statistics relevant to the distribution of three- 
prime Carmichael numbers are given, with particular reference to the conjec- 
ture of Granville and Pomerance in |10| . 



1. Introduction 

A Carmichael number n is defined by the property that n is composite and 
= a (mod n) for all a. For n = HiLiPi"' j with pi prime, Korselt in 1899 [11] 
stated that = 1 for all i and lcm(pi — l,p2 — 1, ■ • • ,Pd — 1) divides (n — 1) is a 
necessary and sufficient condition for n to divide (a" — a) , but he did not exhibit any 
such number n. In 1910 Carmichael f3j showed that the above condition required 
d > 3 and all pi to be odd, and gave four such numbers, the smallest of which was 
561 = 3 • 11 • 17. In 1912 [3] he amplified his remarks and extended his list to fifteen 
such numbers, including one with d = A (although very curiously he reconsidered 
and rejected 561!) 

Korselt's criterion, stated above, is the basis for much of the theory on Carmichael 
numbers and for algorithms to find them, including ours. For a background on 
Carmichael numbers and previous counts of Carmichaels up to increasing upper 
bounds see Ribenboim [16], counts which have now culminated in Richard Pinch's 
up to 10^^ [13]. Our list up to lO^'' for c? = 3 may be found on the website of the 
Cambridge University Department of Pure Mathematics and Mathematical Statis- 
tics [7]. 

In addition to Pi and n already mentioned, Korselt's criterion spawns numerous 
other variables, some of them specific to c? = 3, and various relationships and 
bounds connecting them are of value in constructing algorithms to find Carmichael 
numbers as well as being of interest in themselves. In the next three sections we 
review and extend such relationships and bounds. 

2. Notation; variables; Korselt factorisations, numbers and families 

2.1. Notation, KN's and K-families. Because of the form of the Korselt crite- 
rion, we shall find it convenient consistently and exclusively to use the abbreviation 
x' := X — 1. So we have {xy)' — xy' -\- x' = x'y' + x' + y', etc. We shall also consis- 
tently use the notation which we define during ^ without repeated explanation. 
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Let n = Jli=iP«i where 1 < pi < p2 < . . . < pd and d > 3, be a number for 
which the factors pi satisfy the condition p/ divides n' , and define P := 11^=1^ P«i 
p Pi, q := Pd-i and r := pd, so for d — 3, n ^ pgr. Let Pi := n/pi, so n — piPi 
for 1 < i < d, and similarly write n — qQ — rR. 

Then n' = {piP-i)' = Pi Pi + P/, so n' = (mod p/) gives Pi = (mod p/) for 
1 < « < and conversely. 

Thus there exist integers \i such that Pi = XiPi . Also if Ad = 1, then obviously 
Pd = Pd = Yli=i Pi is necessarily composite (usually we shall assume that Xd> 2). 
We write D := Xd-i , E := Xd , so 

(2.1) Q' = Dq' and R' = Er' 

So far in substance but not in notation we follow Carmichael, if pi are all odd 
primes. But both for algorithms and theoretical results it is often necessary to 
consider sets of numbers n with a factorisation which satisfies the Korselt divisibility 
criterion without all (or any) of the factors necessarily being prime: such a number, 
n as above, together with the particular Korselt factorisation, we shall term a 
Korselt number (abbreviated to KN , or KdN if its Korselt factorisation has d 
factors) if, for all i,pi is odd and E > 2. I have not established whether it is possible 
for a number to be a K^N in more than one way, but the Korselt factorisation will 
always be apparent from the context. Clearly if every pi is a prime then n is a 
Carmichael number, which we shall abbreviate to CN or CdN in like manner. We 
shall also consider certain families of KNs (if-families or ii'd-families) of the type 
"'(^) ~ riiLi Piit): where n and pi are polynomials over the integers and t is a non- 
negative integer parameter. It seems plausible to conjecture (with Schinzel, see 
page 91 of [T^) that any i^T-family will contain an infinite number of CN^s unless, 
speaking loosely, there is some obvious reason why (almost) all members have at 
least one composite pi. 

2.2. Chernick's universal forms. The best known i^T-families are the "universal 
forms" described by Chernick in 1939 [5], and it will be helpful to summarise his 
theory in our notation. Let n be any KN. Then we have 

d d 

(2.2) n' = l[{p./ + l)-l = l[pUj2{pi'p2' ■ ■ ■Pd-i') + . ■ . + Y,Pi'p2' + Y.Pi' 

i=l 1=1 

Let H = gcd[{pi'}] for any particular j. Then from (|2.2p n' = pj' (mod H). But 

since n is a KN, clearly for i ^ j we have H\pi'\n', so n' = (mod H). Hence 
Pj' = (mod H), whence H = gcd [{p/}]. So if Ai := Pi /H, any set of {d—l)Ai's 

l<i<d 

are coprime, i.e. 

(2.3) gcd[{AJ] = 1 , for l<j<d 

Also if L := 1cm [{p/}] and £ := 1cm [{Ai}], clearly L = £H. Combining 

l<i<d l<i<d 

Pi' = A,H, L^£H and n' = (mod L) with (j^ . we get 

H"-^ V A,A2 ■ ■ ■ Ad-i + H"-^ V Ai^2 • • • Ad-2 

(2.4) ^ ^ 

+ --- + i?^AiA2 + ^^i = (modi) 
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Suppose now that we are given any set {Ai, A2, . . . , Ad} satisfying (j2.3p . then con- 
gruence ()2.4|) is always soluble for H when d = 3 (see below, following ()2.5|) ). but 
not necessarily when d > 3. If Ho is any solution, then so is Ht = Ho + t£, so we can 
choose Ho to satisfy 1 < Ho < £■ Then if we take pi = AiHt + 1 = Aitt + AiHo + 1 
and n = OiLi^'ii satisfies the Korselt criterion for all t, and with certain precau- 
tions yields a i^Tt^-family corresponding to each basic solution Ho (precautions: our 
definition of a KN requires (i) all pi are odd, and (ii) E > 2: for (i), if £ is even 
and H is odd, from (12.41) we get 

d d 

^(AiA2 • • • Ad-i)+- ■ A1A2+Y, M = n(^'+i)-n ^^-1 = (^^'i 2) , 

i=l 1=1 

whence since £ \ H^Li ^-i^d £ is even, OiLil^i + 1) = 1 (mod 2) and so for all i, 
Ai = Q (mod 2), contradicting (|2.3p which holds by hypothesis; so if £ is even then 
H is even and all pi are odd; but £ is odd iff all Ai are odd, so Ht is alternately odd 
or even as t increases, and then the iC-family will be generated by the parameter u 
where t — 2u or 2u+l according as Ho is even or odd; while for (ii), i? > 2, it may 
be necessary to exclude t = from the family). These i^- families are Chernick's 
universal forms, of which the best known arises from (^1,^2,^3) = (1,2,3) with 
Ho — 6 and then as above n = (6t + 7)(12t + 13)(18i + 19) (Chernick equivalently 
has (6M + \){12M + 1)(18A/ + 1)). 

For d = 3, let A:^ Ai,B -.^ A2, C A3. Then, for any Ky,N , from (US)) A, B, 
C are pairwise coprime, £ — ABC, and so from (|2.4p there exists a positive integer 
F such that 

(2.5a) H{AB + AC + BC) +A + B + C = FABC, i.e. 

(2.5b) F = + 1 + 1) + + ^ + is a positive integer. 

Also given any pairwise coprime A, B, C, then ^ AB and ABC are coprime, 
so (2.5a) has a unique solution for F and Ho, and we get a K^-iaimly as described 
above. Putting p/ — AiH in (2.2) with d = 3, with (2.5a) we have 

(2.6) n = ABCH{H^ + F) 

Theorem 2.1. For any K3N, I < F < 2H . 

Proof. Based on (2.5b), write F = F{A,B,C,H). We have iJ > 2 and B > 2. 
Obviously F > 1. We consider two cases: (a) B > 3, (b) B = 2. 

19H 2 

(a) For B > '3, F{A, B,C, H) < F(l,3,4,i7) = -j^ + - < 2H ioi H > 2 , 

23H 5 

(b) B=2. For C > 7 we have F(l, 2, C,i/) < F{1,2,7,H) ^ -j^ ^ 7 - 
2H for H >2. Also 

F(l, 2, 5, H) < F(l, 2, 3, H) = + 1 < 2i/ for > 6. 

6 

But for {A,B,C) = (1,2,3) or (1,2,5) we have Hq — 6, whence the result, with 
F = 2H only for n = 7 • 13 • 19 = 1729. □ 

Note: {A,B,C,H) = (1,2,7,2) yields n = 3 • 5 • 15 = 225, not a K^^N since 
E^l. 
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2.3. The equation(s) of Beeger, Duparc and Pinch. For any KN we have 

n = Pgr = qQ = rR, and hence Q ~ Pr and R ~ Pq. 
Then from Dq' ^Q' ^ {Pr)' ^ Pr' + P', so 

(2.7) Dq' -Pr' = P' , and similarly Er' ~ Pq' = P '. 
Solving simultaneously for q' and r', and writing 

(2.8) A:=DE-P^, we get 

(2.9a) q' = P'^P^^) and 

(2.9b) ,^P'{P + D) 



A 

Beeger for d = 3 in 1950 [2] and Duparc for d > 3 in 1952 [9] gave (EJa]), and 
Pinch fT2^ bases his first algorithm on (I2J1) . From ((2J1) . A > 1. Also Duparc 
showed 

Theorem 2.2. For any KN, 2 < E < P - I. 

Proof. From the definition of a KN, E >2. Also r — q — 1 > 1. Hence from ()2.7|) 

Pq' + P' _ Pr' - Pjr' -q')+P' _ P(r - g - 1) + 1 ^ 

whence the result. □ 

The following equations based on (|2.8p and (|2.9p will also be useful in §3 and §4. 
Define s := P — E and r/ := D — P — s, so 

(2.10) E = P - s and D = P + s + r] 

Then from Theorem [221 1 < s < P - 2, and from we have 



A + A + 

(2.11) A = rj{P — s) — , whence rj = — — — — — and so 77 > 1 

P — s E 

Hence 



(2.12) s'^ + r]s ^ TjP - A, so A < r/P, and s = W yyP - A + ^ - ^ 
Hence if 



(2.13a) e := ^4(77P - A) + r;2 , then 

(2.13b) * = ^ = ^"^' D = P + 

Then (1^ becomes 



(2.14a) q=£l(2P+l—)+l , 



A V 2 

P'/ ^ ri + e 



(2.14b) ^ = ^(2^+^)+^' 

P'2 f / „x2 5)2^ p/ 

whence = — |(2P+^j _ _ j + _(4P + r;) + 1 
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and using ()2.13ap for 9"^ we get 

, p/ 2 pt . 

(2.15) n = Pqr = p| — (4p2 + ?;P + A) + — (4P + ??) + l| 

Subject to the solubility of certain congruences, a particular choice of A and rj 
leads to one or more i^s-familics via p.l2l[2T3ll2.14|) . using ()2.12|) with s as a first 
parameter. 

2.4. Relations connecting K^^N variables. These will be needed in §4. 

R' (pq)' p'q'+p' + q' ABH + A + P 
We have E = — ^ = = , 

since pi' — AiH , and hence 

(2.16) ^ ^ ABH + A + B ^ Bp + A 

E E 

Combining this with ((23a| we get ABCF = CE + CH{A + B) + C, i.e. 

(2.17) ABE = {A + B)H + E +1, whence, writing 

(2.18) G:=AF- H, we have 

AH+E+1 p+E 



(2.19) B 



G G 



.^^^ , Pip + E) AHlp + E) Aip + E) 

Also using (Eli). BH = q' = ^ J ' = \ ' , so B = J ' , 

t, M V, A A A 

whence from (I2.19p . 

(2.20) A = . 

So G > 1, we have A > 1, and we now show 

Theorem 2.3. For any K3N, I < G < 2H and I < A < 2p - 2. 

Proof. We have p < q' and, from Theorem 12. 2[ E < p' . Hence, using (|2.20|) 
and ((2Ja|) . 

AG^A^P-^^P±^<^-^^^^^ = 2p-Z+^-^2AH-l+^-<2AH = 2p' , 
q' p p p 

from which both statements in Theorem 12.31 follow . □ 



We observe that (|2.181 12.191 12.16P enable us to express in turn A, p, P, q, C, r 
and n in terms of E, F, G and H, which we shall exploit later. 
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2.5. Bounds, variables and challenges. The next two sections are concerned 
witli finding inequalities y < f{x) showing upper bounds for y given a;, where x 
and y are variables connected with KN^s and hence CiV's. Here / is an increasing 
function and, usually, x is P, p or n: for example, we shall show that an upper 
bound for r given n is given by r < . Obviously, if exists, any such 

relation is equivalent to a lower bound of f~^{y) for x given y. Also, if a; is n and 
we are looking for all CN^s less than some large X, we have y < f{X). Our symbols 
have been defined as integer variables connected with KN^s, but sometimes we shall 
treat them as real variables obeying the same inequalities as the integer variables. 
If possible we shall exhibit if-families for which the bound is attained, or failing 
that some CiV's or KN^s for which it is approached. Some of these bounds were 
used in executing our algorithm for C^N^s (see §5), although invariably a weaker 
(and more easily proved) bound would have sufficed. 

Challenges: occasionally I offer a challenge to find a C^N or K^N satisfying 
certain conditions, and I would be very interested to receive successful responses at 
my address at the end of this paper. 

3. Bounds for KN variables with d > 3 
3.1. Upper bounds given P for q, r and n. 
Theorem 3.1. (Duparc's theorem) For any KN, 



,<(p-l)(2P+i-^pJ)+l. 



Proof. From (|2.9ap and Theorem 12.21 if A > 2 we have 

(3.1) ^,^P'iP + E)^P'iP + P')^ . I 
^ ' ^ A - 2 V 2 
Fix A and P, with A < P, and regard 6 — 9{ri) and q — q{if) as functions of a 
real variable 77, defined by (|2.13ap and (|2.14ap respectively, with 77 — rjo for the KN 
under consideration. Then we have 

(3.2) U^~0) = l- , "-"'"^ <0, 

V(?7 + 2P)2-4(P2 + A) 

so as 77 increases, r\ — Q and hence q both decrease, and hence q{r\o) < q{l). We shall 
use (13?2]) for general A later, but with A 1, q{l) = P'{2P+\-J P - |) + 1. But 



for P > 3, P'(P- |) < P'{2P+\ - ^P - |), so with the result follows. □ 

We note that for this maximal q situation, with A = 77 = 1, (|2.15[) gives 
(3.3) 71 = A^2(P) :=4P'^-7P'* + 7P^-4P2 + P 

Also from (|TT^ [TT^ I^Ti)) we get 
(3.4) 

P=P2(s) :=s2 + s + l , 6' = 2.s + l 



(7 = g2(P) := ( P - 1 ) ( 2P + i - Jp - ^) + 1 =92(5) 25* + Ss^ + Ss^ + 2s + 1 



= P2(P) :=(P-l)(2P+i + A/p-^)+l =r2(s) := + Ss^ + + 3.s 
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and we have the q-maximal A'3-faniily n — 712(5) := P2(s) ■Q2{s) ■r2{s). A computer 
search by Ian WiUiams (see ^J7]) for 1 < s < 4906 found just 12 g-maxuxial CsN's, 
with s = 1, 2, 6, 90, ... up to 3654, and only one g-maximal CdN with d > 3, namely 
the CqN n2(1493) = 7 • 19 • 31 • 541 • 9947309489407-9953972118361 ^ 2.2 x 10^2, 
with P = 2230543. 

Theorem 3.2. For any KN, r < RsiP) := ^{P - 1)(P + 1)^ + 1, with equality iff 
q^QsiP) :=p2 + p_i. 

Proof From dm and (H^bJ, 

A A V E J \A AE EJ ' 

But A > 1 and > 2, so 

/<P'(p+^ + l)=ip'(P + lf , 

with equality iff A = 1 and E = 2, which from (|2.9ap and p.ip occurs iff 
g' = P'(P + 2), whence the result. □ 

Theorem 3.3. For any KN, n < N^iP) 1{P^ + 2P^ - P'^ -P^ + 2P'^- P), 
with equality iff q ^ Qz{P) = P'^ + P - 1 and r = RsiP) = j{P - 1)(P+ 1)^ + 1. 

Proof. For A > 2, from g < _ |P + | < Qz{P) for P > 3, and from 

Theorem [321 r < R^{P), so n < P ■ Q^{P) ■ R^{P). 

For A = 1 and given P, with 77 and n real variables, from (|2.15p 71 is greatest 
when 77 is greatest. But from (|2.1ip 

A + + 1 



and 1 < s < P — 2, giving maximum rj = ^(P^ — 4P + 5) when s ~ P— 2, i.e. E ^ 2 
from ^JQ. But from Theorem E21 for A = 1 and S = 2 we have q = Qj.{P) and 
r — R^{P), so n < P ■ Q^{P) ■ R^{P)^ which multiplies out to give the result. □ 

Beeger [21 for d = 3 and Duparc 9J for any CN proved results similar to Theorem 
13.11 and weaker than Theorem 13.21 They were chiefly concerned to show that the 
number of CN's for given P is finite. Swift stated the first result of Theorem 13.31 
for d = 3 in 1975, but his proof is not pubHshed in [T7] . 

For a ii'3-family attaining these upper bounds for r and n given P, we simply 
put P = 2i + 1 in 77 = N3{P) = P ■ QaiP) ■ RsiP). We note that Q2(3) = 
(53(3) = 11 and ^2(3) = ^3(3) = 17, so the smallest CN (and KN, easily shown), 
77 — 561, uniquely is g-maximal, r-maximal and 77-maximal. N^iP) gives CsN^s for 
P = 3, 5, 31, 41, 83, 131, . . .; and iV3(65) = 5 • 13 • 4289 • 139393 = 38860677505 is the 
smallest r-77-maximal C4N for given P. A computer search over 3 < P < 132425 
found 178 C3iV's, 18 C4A^'s, 29 C^N's and 9 CqN's which are r-77-maximal, the 
largest of which is the C3A^ with P = 131711. 
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3.2. Upper bounds given n for P, q, r. For an upper bound for P given n, we 
do not attempt to improve on the obvious: 

Theorem 3.4. For any K^N, P < rT^ 

Proof. We have pi < p 2 < ■ ■ • < P rf-2 < q < r (with obvious modification for 

d-2 

Pi < P^l9, SO pd-2 > P^ and P = — < ^5— < —. 

,^1 Pd~2 P^ 

Hence < n and the result follows. □ 

It seems plausible that an upper bound for q given n should correspond to the 
Theorem 13. II bound for q given P, so in terms of l\3.'S\i and p.4p we should have 

Theorem 3.5. For any KN, q < Q2{N2^{n)) = (72 ("-2^ ^(«)), with equality iff 



Q,[P) = [P-l][2P+\-^P-\ 




Explicitly, 



Proof. liV ^ \n), the substitution V ^ + v + \ with the algebra of (|3^ 
and p.4p establishes the equivalence of the functions (52^2^^ and (72 "-2^^, and if 
A = 77 = 1 it is immediate that q = Q2{N2^{n)); also A = 77 = 1 for 71 = 561, 
the only KN with P = 3. So we assume KN^s with A?7 > 2 and P > 5, and we 
consider two cases: (i) A > P', (u) 1 < A < P - 2. We write A := 

(i) A < 1. From (|2.13bp and (|2.14ap . q' = A(2P-s) < 2AP, and then from ([^1^ 
n = P{A2(4P2 + ?7P + A) + A(4P + ry) + 1} > 4A2p3 > 4A2(|^)3 = ^. Hence 
g'"^ < 2An < 2n, so q < (2n)3 + 1. Let v := n2^{n), so n — n2(v), and it is easily 
verified that (2n2(w))3 + 1 < q2{v) for v > 1.1. But 712(1. 1) < 963 and every KN 
with P > 5 has n > 1105 = 5 • 13 • 17 and hence v > 1.1, so q < q2{n2^{n)) as 
required. 

(ii) A > 1. Let Ug be any particular KN with associated ii'- variable values 
Po, <7o, "-o, Ao, 770, 610, So and Ao = > 1. We define s, 9, q, r and 77 as functions of 
independent real variables P, A,77 by the formulae of (|2. 12112. 15]) over the domain 
3 < P < Po, 1 < 7; < 770 and 1 < A < Ao with A < P (ensuring 9 e 5R), and 
we write n — n{P,A,ri) with Uo = 7i(Po, Ao, 770), etc. Then, if 770 > 1, over the 
interval 1 < 77 < 770 from (|2.14ap . ((3?2)) and ((2?T5)) we have || < and |^ > 0, 
whence qa ■— g(Po, Ao,l) > qo and ria := 77(Po,Ao,l) < 77o, with equality iff 
rjo = I . In what follows we use the suffices "a" and "b" to correspond to "states" 
(A,7y) = (Ao, 1) and (A, 77) = (1,1) respectively. 

Then, if Ao > 1, for 1 < A < Ao we define the function P P*(A) imphcitly 
via the relation 77(P, A, 1) — Ua, so ria = ni, and 

(3.5) p|^(4p2+P) + 5^^ + l| = P{A2p(4P+l) + 5AP+l} = 71a, constant, 

and as A decreases from Ao to 1 we see from the first of these expressions that P 
steadily decreases, whence from the second A steadily increases. So, for Ao > A > 1, 
= A > Ao > 1, whence P - A > 1, so 6* = ^4(P - A) + 1 > ^5, ensuring that, 
via (I2l3)l and (jTTi)) with 77 = 1, 9*{A) := 6'(P*(A), A, 1) = 9, and, similarly, 
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q*{A) = q and r*(A) = r are defined for 1 < A < A^; also F > A + 1 > 2, 
so certainly N2{P) is an increasing function over tlie relevant domain. Now, for 
> Ij write Pf, := P*{1) ~ N2^ina), so since by hypothesis Po > 5 and Aq > 1, 
we have 

N2{Pb) ^n{Pb, 1, 1) - na - n(Po, A„, 1) 

=Po{\l{4P^ + Po) + 5\oPo + 1} > + + p^> 655 > 561 = iV2(3), 

whence Ph > 3; but P > P*(l) = P;,, so P > 3. Thus from the above we have 

(3.6) For 1 < A < Ao , A > 1 , P > 3 and 6l > VS. 

Also let Sb ■.^P2^{Pb) = 5(V4P6-3 - 1) and qt := Q2{Pb) = q2{sb), so if A^ = 1 
then qb = qa- 

Next we shall show that, with r] — 1, ^ ^ 1 < ^ < ^o, and then 
since (32,92,^2 and ri2 are all increasing functions over the relevant domains, we 
shall have q*{Ao) < and so qo < qa = g(Po,Ao,l) = q*{Ao) < q*{l) = 

q{Pb, 1, l)^qb^ Q2{Pb) - Q2{N^\na)) < Q2{N^Hno)), and likewise go < % = 
g2(s6) < 92(n-2^^(?io)), with equality throughout iff Ao = rjo = 1, i.e. iff = Q2{Po) 
as required. 

From the first of the two expressions for Ua in (|3.5p we have 



dA A2^ ' ' ^ 'UA A3 

5(2P-l)dP 5PP' 



that is 



A dA A2 
(3.8) 

d(gr) _ rP'(16P2-5P-l) 5(2P-l)-,dP 2PP'2(4P+1) 5PP' 

dA ~ I A2 ^ A J dA A3 A?~ 

dP d(qr) 

But Pqr — Ha, constant, so qr — - + P — = ; 

dA dA 

substituting from (|3.7[) and (|3.8|) gives 

fP'2p(4P+l) 5PP' PP'(16P2 - 5P- 1) 5P(2P-l)-idP 

I A2 ^ A2 + A JdA 

_2p2p'2(4P+ 1) 5p2p/ 



A3 A2 

A2-1 dP 



whence 



(3.9) a|2P'(10P^-4P-1) + 5A(3P-2) + — |— = 2PP' 2(4P+ 1) + 5APP' 



Also with 7? = 1, (|2.14ap gives 2q = ^(4P + 1 - 6*) + 2 and from (|2.13ap 
6»2 = 4P + 1 - 4A, so 6* = i(4P + 1 - 4A), 

de 2/dP \ , ^dq fldP P'\/^ ^ \ P'i.dP dO \ 
dA = ^(dA-0 2di = (AdA-A^)(^^+^-^)+A(^dA-dA) 



2 dq ^^dP ^,\(,^,, AP +1 ^ AA\ ^ ^ _^,( ^dP 2dP^2 



Hence 2A^-^ = (A--P'j(4P+l-^+-j+AP'(4---- + - 

, ,,dg ,f 6P-1 4A -I dP ,/ \/ 1\ 2AP' 

(3.10) 2A^^ . a{8P- 3 - ^ + - ^'(4^+ - ^) - ^ 
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Writing M := 2P'{IQP'^ - 4P - 1) + 5A(3P - 2) + ^, from and (IXTUl) 



p 

we get 

-2A'M^ = {P'(4P+ l) (l - ^) + ^^}{2P'(10P2 - 4P - 1) + 5A(3P - 2) + ^} 

_|8P - 3 - + ^}|2PP'2(4P + 1) + 5APP'} 

A 

The coefficient of — in this expression is 


P'{4P'(10P2 - 4P - 1) - 5(4P + 1)(3P - 2) + 5P(6P - 1) - 8PP'{4:P + 1)} 
= P'{4P'(2P2 - 6P - 1) - 5(6P2 - 4P - 2)} = 2P' ^{2{2P^ - 6p - 1) - 5(3P + 1)} 
= 2P' 2(4P2 - 27P - 7) = 2P' 2(4P + 1)(P - 7) 

A2 . 

and of — — is 
6 

p'l 4^'+! ^^^^ _ P'(10P2 - 24P- 1) 



^'1 — + 10(3P - 2) - 20P| 



.2A2M--^ = 2P'2(4P+l)(2p2-P-l) + 5AP'(4p2-2P-2) ■ ^ ' 



dA V ' ^ ' ^ / ■ p 

2P^^(4P+1) ^^^, ^ 2AP-^(4P+l)(P-7) 

A2P'(10P2 - 24P - 1) 2A3p' 
H H 

ep 0P 

Factorising terms and dividing by 2P' ^(4P + 1) gives 

A^M ^^^^2P+i . 5^(2P+1) , A2 (4P+1) 



(4P + l)P'3dA ' P'(4P + 1) 2PP'2 

A(P-7) A2(10P2 - 24P- 1) A 



eP' 20PP'2(4P+1) 6iPP'2(4P+l) 

From (13SD P > 3 and 6* > \/5, so M > and also lOP^ - 24P - 1 > 0, whence 

W>2P+l-'-^+'-^^^ + ^ 
e A(4P+1) 

_ (26* - 4)P + 61 - 1 (4P + 1)(P - 7) + 56i(2P + 1) 

^ e ^ A6'(4P+ 1) 

(4P + 1)(P - 7) + 5%/5(2P + 1) _ 4P^ + (lOVS - 27)P + 5%/5 - 7 

^ A6i(4P + 1) ^ A6i(4P + 1) ^ 

since the discriminant of the numerator is negative. 

dq (4P + 1)A'^P 'W^ 
Hence — — = < as required, completing the proof of the 

first part of the theorem. 

We now outline a method of expressing (52 (-^^2 as a power series in 

Let [/ := , and ((3?2l) becomes 

(3.11) [/5^p5_Ip4^Zp3_p2^1p 

^ ^ 4 4 4 
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and if wc now put P = U + Bo + 'jf + §i + ---, substitute into (|3.1ip and equate 
coefRcients of U^^'' to find successively Bo, i?i, B2, . . ., we obtain 



(3.12) 
To find 



7 



21 



1 



2787 



20 200U 250[/2 160000C/3 



+ 



P 



and with u ^ VU = ^{/f , from ([XT^ we have 



P 



21 



7 



20 200^2 



1 



P-l 



2 2 
It — - 



250ii4 
21 



P' 



2PP' 



13 



21 



P\ P 



1 



20 200^2 250u4 



and 



1 



5 200u2 250u'i 



21 



1 



5w2 200u4 250^6 
1 29 1 



53 



477 



after applying the binomial series and simplifying. 
Substituting back, we get 

^ „ 4 c, u2 17u 1 19 

Q2{P) = 2u^-u^ \ \ \ ^ ^ + . . . 

' 10 20 5 400w 2000u2 SOOOu^ 

If now the variables all belong to a KN, then q — Q2{P) is an integer and 
q = [2u4 -u^ - 0.1^2 + 0.85uJ , and Theorem [33] follows. □ 

Much more simply, we now establish an upper bound for r given n; first we prove 

Theorem 3.6. For any KN, 



Pi 



+ 4A, 
2A,; 



Proof. Using the notation of ^2.1[ we have 

PiPj ' _ p- 

^ a" ^ ~ 



PiPi 



A,; 



A,; 



whence 



(3.13) 



n 

a" 



Pi 



and XiPi 2 — A,; 'pi — n = 0, giving the first equality. 

Also A, > > 2, so < 1 - ^ < 1, and from ((XT^ we thus have 



(pi- 1)2 < <p,2^ 
Al- 



so Pi 



completing the theorem. 
Theorem 3.7. For any KN, 



□ 



r < 



y/8n +1 + 1 



and r < 



with equality iff E 
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Proof. This follows easily from Theorem 13.61 

n > 561 and n is odd, so with i = d in Theorem 13.61 and Xd = E > S, we have 



In V8n + 1 + 1 
and \ — < 



/n _ n 
£ + ^<V2< V2 V2 - 4 

in conjunction with ii' = 2 in Theorem 13.61 Theorem 13.71 follows . □ 

The smallest CN which is r-maximal for given n but not for given P (see 
Theorem [32]) is 8911 = 7 • 19 • 67. Another of this type is 949803513811921 = 
17 • 31 • 191 • 433 • 21792241, which Pinch in T2] says contains the largest prime 
factor among CN's < 10^^. 

4. Bounds for K^N variables 

4.1. Upper bounds given p for A, B, C. To establish an upper bound for A 
given p, we need one for A given H: 



Theorem 4.1. For any K3N, A < 3iJ - f . 

Proof. Suppose for some K3N that A > \H for some A > 0. Then (|2.5bp yields 
I I \ I 1 1 



F < H 



\XH XH + 1 XH + 2J XH{XH + 1) XH{XH + 2) {XH + l){XH + 2) 
1/1 1 \ /l 2 \ 3XH + 3 



X \X X{XH + l)J \X X{XH + 2)J XH{XH + 1){XH + 2) 

3 iXHlH -l)+AH -i 3 . „ 

'— < - since H >2 

X XH{XH + 1){XH + 2) X 

But F > 1, so putting A = 3 we get A < 3H. Also if A = | then F < 2, whence 
if 3ff < <- 3^^ then F =1. We define a UgA-K^N to be a K3N with A>^, 
and likewise a higA-C^N . All other K^N's obviously obey Theorem 14. II 

So we now put A = 3H — a, with 1 < a < and using F = 1 we shall 
show that for given a, H < 2a} , yielding Theorem 14.11 We write a := —a, 
v4 = 3i7 + a, B = 3i7 + /3, C = 3i/ + 7, cr := r := so ^7 = - r^, 

and S* := ^ a = 2(7 — a. Then for a K3N we have —a — a<P = (7 — t<j, 
so < T < a + cr, and from (|2.5p 



+ _ 27g3 + a)H^ + Q/3 + 9)g + ^ a m 
ABC ^ 27i/3 + 9Ea)i/2 + 3(^Q,^)ij- + Q,^^ ^ ^ISC 

where 

(4.1a) m := TO(iJ) 3(^ + (2 ^ a/3 - 9)i7 + a/37 "51" 

(4.1b) = 3SH^ + (2(7^ - 2r2 ~ iaa - 9)i/ - a{a^ - - 1) - 2cr 

So for any F = 1 iff m = 0, and A = 3H - a > 0, so > |; also for any 

higA-K^N , 7 > (since otherwise 

H H 



1 H 
>- and so F > V 7777 > 1 , but F = 1). 



C 3i/ + 7- 3 ^3iJ + a 

We now regard m(H) and F(iJ) := 1 — —— tttttt . 7 as functions 

^ ^ ^ ^ (3i/ + a)(3i/ + /3)(3i7 + 7) 
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of an unrestricted real variable i7, where a, /3, 7 are real, a < < 7 and a < (3 < ^ . 
Essentially by considering the graph of F{H), we show that m{H) = has a root 
ff* > I = -| iff 5 > 0. We have 

9 

(1 or (X 

m(--) = (a + /3 + 7)y - {2a(/3 + 7) + 2/^7 - 9}- + a/37 - (a + /3 + 7) 

= |{a' - a(/3 + 7) + /37} + 2a - /3 - 7 

= |(a - /3)(a - 7) + (a - /?) + (a - 7) , 

with similar results for to(— 'I) and m(— ■^), whence r7i(— ^) < and m(— ^) > 0. 

First we show that in all cases m{H) = has a root h* satisfying — ^ < h* < — ^. 

^ a 
(4.2) As H^-^- and as -- + , F{H) ^ 00 

Then if /? = 7, as ^ + , F(i7) ^ 00; as 77 ^ -f-, F{H) -00; and 
F(iJ) is continuous over the interval so there exists h* G (^3 ; ~f ) with 

F{h*) — 1 and m{h*) = as required; and similarly if /? = a. Also if a < /3 < 7 and 
(a) m(-|) ^ 0, then as H ^ and as F ^ -f-, F{H) -cx), and F{H) 

changes sign as H increases through the singularity at H = — ^, so again there 
exists h* with F{h*) = 1 and m{h*) = in at least one of the intervals {~-^, — f ) 
and (-|, -f ), and so in (-^, -f); while (b) if m(-|) = then h* = -|. Hence 
m{H) has at least one zero h* £ (— ^, — as required. Then, if = 0, m{H) is 
linear and h* is the only root. But, if S* > 0, m{H) is quadratic and as H 00, 
F{H) 1-, so with dig) this gives H* > -f such that F{H*) = 1 and thus a 
unique second root H* of m{H) — with H* > |. Similarly for 5 < 0, the second 
root H* satisfies H* < — ^ < 0, and hence iJ* > | iff 5* > 0, as required (in fact 
H* > ^, else F{H*) > 1, as is easily seen). 

Assuming henceforth that S > 0, from the above argument for iJ > | we have 
F{H) > F{H*) = 1 iff iJ < H*. So if, with obvious notation, for (aj,/3j,7j), 
ai = ^2 = -a, and for every iJ > | we have Fi{H) > F2{H), then Fi(iJ|) > 
^2(^2) = 1 = Fi{Hi), so H2 < HI. In particular we deduce that 

(i) if ai = a2, /3i < /32, 7i < 72, with at least one strict inequality, then 
H'2 < HI; 

(ii) if ai = a2 = —a and Si = S2 — S (so a\ = a2 — (y — ^^)i but ri > r2, 
then 7?2 < -H^ii i-^ for fixed a, and u, i/* increases as r increases. This follows 
because (a, /3, 7) = (a, ct — r, cr + t) and 

n^)^'^Ei + Ei,-f + ife{(H + 3)(^ + ^) + '}^ 

but S + C = 6i7 + 2ct and = (3i? + cr)^ - r^; so for > f , Fi(i?) > F2{H) 
and iJj* < 

We write H^{a,S,T) i?*(a,/3,7) iJ*, and /i(/3,7) := iJ*(-l,^,7). We 
observe that for a \ngA-K^N , H ~ H* , and that A, B, C pairwise coprime and 
H even requires that no two of a, /3, 7 are even; so, since f3 + ^ = S + a, S odd 
requires a, f3, 7 all odd. We now show that for any higA-K^N, H < 2a^, considering 
cases (a) a = 1, (b) a > 2, 5 = 1, and (c) a > 2, 5 > 2. 

(a) Any K3N with a = 1 is a higA-K^N. So a = -1, and for /3 = 0, since 5 > 
and B, C are coprime, 7 > 5 and from (j4.1a|) and (i) above we have 
/i(0,7) < /i(0, 5) ^ 1.68 <2<H. Also for /3 > 2, by (i) we have 
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h{l3, 7) < h{l, 7) < h{l, 3) - 1.63 < 2, covering all cases except h{l, 2) - 2.14 ^ 2N. 
Thus there are no K^N's with a = 1. 

(b) We have 5 = 1, odd, so a is odd and a > 3. Then maximum r for given a 
occurs when {a, (3, 7) = (—a, —a + 2, 2a — 1), giving r = ; and from (ii) above 
H < g ■— g{a) := H^a, 1, ^^^^), which from (|4.1ap is given by 

3g^ - (6a^ - 8a + 13)g + 2a^ - 5a^ + 2a - 1 = 0, whence 

g = i(6a^ - 8a + 13 + \/i6a^ - ISOa^ + 280a2 - 244a + 181), 
i.e. 6 

(4.3) H < g{a) = ^{6a^ - 8a + 13 + ^/(6a^ - 10a + 15)2 + 68a - 44) 
< 2a2 for a > 3. 

(c) First we find H^{a, S,a + a) = H*{-a, -a, 2a + S), which from (|4lb)) is the 
root H* of 35i?2 „ (20^ + 8acr + 9)H + + 2a'^a -2a + a = 0, which has 

discriminant = {20^ + Sacr + 9)^ - 12(2cr - a)(a^ + 2a'^a - 2a + a) 

= 16a^ + 32aV + Wa^a^ + iSa^ + 96aa + 48a^ + 81 

= 16(a2 + 3)(a + cr)2 + 81 and hence 

1 



(4.4a) i? t (a, 5, a + cr) = (2a2 + 8acr + 9 + Jl6{a^ + 3){a + a)^ + 81) 

65 

(4.4b) = ^{2a^ + Ma + 5) + 9 + ^y4{a^ + 3)(3a + 5)2 + 81) 

65 



. ^ a2 2a 3 // , \fa 1\2 9 

(^•4^) = 5+y+25 + Vr+')(5 + 3) +45^ 

So from (i) or (|4.4cp for given a, H^a,S,a + a) decreases as 5 increases, and 
hence for any big^-i^siV with 5 > 2 and a > 2, also using (ii) we have 

H = H\a, 5, r) < H\a, S,a + a)<H^ (a, 2, ^^^) , 
and from (j4.4bp we have 



i7t(a,2, 



3a + 2\ ^2 f oA/a 9 3 2 2a 3 

< 2a2 W a2 + 3 - + - + — < -a2 - — - - 



2 J V V y V2 3/ 16 2 3 4 

a-* a^ /3 1\ , 1 9 4 „ . /4 9\ , 

^- + - + (- + -)a- + a+-<-a^-2a-V(---)a2 + a 

<^ 6a'' - 7a^ - 8a2 > 1 

^ a^(4a - 7) + 2a2(a2 - 4) > 1, which is true for a > 2. 

Thus for aU bigA-A'giV's H <2a'^, and Theorem O follows. □ 

It is easily shown that H < g{a) < 2a^ in case (c) as well as case (b) and, for 
H>6, 5-i(iJ) > ; hence for all ifgiV's 

(4.5) A<3H^g~\H), 



which for > 6 is a slightly stronger but less convenient result than Theorem 14. II 
It is also not onerous to extend the approach of case (a) for finding all K^N^s with 
a = 1: for a = 2, again there are none (although 9801 = 9 • 11 • 99 obeys the Korselt 
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divisibility criterion); but for a = 3 there are two, the higA-CsN 7 • 23 • 41 = 6601, 
and the higA-KsN 

(4.6) = 547 • 575 • 659 = 207271975. 

Actually from (gS]), .g(3) = 14, giving rise to n* with (A,S,C) = (39,41,47) and 
{a, 13, 7) = (-3, -1, 5); and 39 = A < 3iJ - ^ 39.35. n* is the only KgiV with 
equality in (|4.5p . since g{b) is irrational, and for a > 7 we have 



6a^ - 10a + 15 < V(6a^ - 10a + 15)^ + 68a - 44 < 6a^ - 10a + 16, 

so from (|4.3p g{a) is irrational. Gordon Davies (see ! j5.ip did a computer search for 
big^-XsA^'s with S = 1, using (a,/3,7) = (— a, — a + 2i, 2a — 2t + 1), for odd a up to 
1239 and 1 < i < no more were found, and H* was rational only for a — 151, 
t = 89, giving H* = 13067i. 
In like manner to (14.41) we find 



H\a, S, 0) = ^{'iaa - 2a^ + 9 + ^4(0-2 + 3)(a + a)2 + 81 ) 
6d 

1 



= T77^{{a + S)i3a - 5) + 18 + ^[{a + Sf + 12](3a + Sf + 324}, 
SO for fixed 5' and large a we havei?t(a,S',0) - andi7t(a^5'^ 

a+cj) ^ ^ , whence 

(1 + o(l))f^ < i?* < (1 + o(l))^. Hence with A ^ and A = 3H - VXH, 
foi a S we expect ^ < A < 2S'; we may regard ^ < 3 and A > ^ as different 
measures of closeness to the bound of Theorem 14.11 For n < 10^^, there are only 
71 bigA-CaA's, and just 11 with ^ > 2; the two largest ^ values are about 2.342 
and 2.683, and only these two bigA-Cs^V's have a > S. With ^ - 2.683, we 
have nt := 835327 • 893359 • 1117117 = 833645090806507981, with {A,B,C,H) = 
(1497,1601,2002,558), so a = 177, 5 = 78 and A - 56.15. For n*, from ((TO . 

— 2.786 and A = I was able to find marginal improvements on the bound of 
Theorem 14. II for sufficiently large a (e.g. (|4.7p below), but none implying A > ^ for 
some > ^■ 

Challenges 1(a): Find a higA-KsN with a > 3 and S ^ 1. 

1(b): Find a higA-CsN with n > 6601 and A < 10. 



Theorem 4.2. For any K3N, A < ^3(p - 1) - ^y^- 

Proof. Using the notation and from the above discussion of Theorem l4.11 Theorem 
holds for n* and for n = 6601, and hence for a < 3. For S > 2 and a > 4, from (|4.4bp 
we have 

3a + 2 , 1 



H < H^{a,2,— — ) =— (6a2 + 8a + 9 + ^/4:{a^ + 3)(3a + 2)2 + 8I) 



1 



< — (6a2 + 8a + 9 + ^A{a + l)2(3a + 2)2) 

2 3 13 „/ 1\2 . 

=a + 2" ^ 12 ^ V ^ 2/ ^"^'^'^ a > 4 
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and for = 1, a is odd, and for a > 5 from (j4.3|) . H < g{a) < 2{a — |)^. Hence for 



a > 4, for any H < 2{a - i)^, so a > W -j + 5 and 



(4.7) A<3i?~y|-i 

Sd' I v' 1 v' 

Let f(x) -.^ J— - X - \ . Then since H = ^, from gJl) f(A) > 0. Also 

X y 2x 2 A 



/'(.)=-|:-l + y|3^-^(6v/V-V^)-l<0 forO<.<p', 

so f{x) decreases as x increases over this interval, which certainly contains A. 
Writing p := -^3p' and p, :— ^ — we have 



^ V6a* 2 

24 



pM(p-#)'-A^n 1 p^ 1 



2 2 

But p > 3, so p > 1.565, p > 1.459, ^ < 1.15, — < 0.00168 

and f{fi^) < 0, whence since f{A) > and /(a;) decreases with x, A < p?, giving 
Theorem □ 



For n* (see ^^), A = 39 < ^/W - \\Jt2 - 39.1736. 

In seeking an upper bound for B given p, we know from Theorem l3 . 1 1 that we can 



have q' = p'{2p — — | + 5), with H = p' and B = 2p ~ yp — | + i; we show 
that it is posssible for B marginally to exceed this: 



Theorem 4.3. (a) For any K3N, B <2p- ^p- \ + 

(b) For any C3N, B < 2p - .^^^1 + ^ 



Proof. Write S=yp— j — ^. We seek B > 2p — S, and we consider two cases: 

(i) G > 2. Then from Theorem [121 < p - 1, and from 

(ii) G = 1. Then from (j^T^ S = p + and from (PTU)) E = p-s, so 

(4.8) B = 2p- s 
Also from A = AG = A, whence from ([2TT|) 

(4.9) , = ^ 

p — s 
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and using ([2Jb| and ^JU^, CH = r' = ^^^P'^ ^ go 

(4.10) C = p + D = 2p + s + ri 
Since H is even, from ()2.18|) 

(4.11) H ^ AF - 1 and A and F are both odd, and also 

(4.12) p = AH +1 = FA^ - A + 1. 

We now consider the two sub-cases (a) 77 > 2 and (b) rj = 1. 
(a) If 77 > 2, from we have s"^ + A > 2{p ~ s), whence s'^ + 2s > 2p - A 
and s > ^2p - A + 1 - 1 > V'2p^~7^^^n - 1, using Theorem [T^ 



But \j2p— VSp' + 1 — 1 > S" reduces top+|>yp— | + \/3p^, which is true for 

p > 3, so s > 5 and from ((48)l B <2p~ S. 

(b) If 77 = 1, from (|49)) + A = p - s, whence from ((4T2)) 

(4.13a) = s(s + 1) = FA^ - 2A + 1 = F'A^ + and 



(4.13b) ,^^FA^-2A-,\-\ 

From (|4.13ap we note that F — \ imphes A' ^ = s(.s + 1), which is impossible. So 
by (|4.1ip for odd F > 3, remembering that a; = + 1 and using (|4.12p again, we 
have 



FA^ -A + --\ FA^ ~2A+- = 



^ V 4 ^J{F'A^ + {A- i)2 + .^F'A^ + A^ + \ 
Hence if = 1, s = as for Theorem |3lll but if A > 3 (A is odd by (SU)) then 
S-s= , > 



F'{i + + (1 + + v^'(i + i^)^ + 1 + ii- 

From this we deduce that < S — s < —K= < since F > 3, 
whence 2p— S<B = 2p— s<2p— S + giving Theorem 14. 3r a'). 

Also we see that for given F, if there is an infinite sequence of K^N^s with 

increasing A values, then 2p — S -\ 7^ — _B^O+asA— >cx). There is such a 

sequence iff (|4.13aP has an infinite number of solutions for s and A which result in 
pairwise coprime A, B and C, and E > 2. Such solutions we call acceptable. With 

(4.14) ^ 2FA -2^2H, 9 = 2s + l, (|4J3a|) implies 

(4.15) (j)^ ~ F9^ = A - 5F, 

and from the theory of quadratic forms a necessary condition for this to have a 
solution is that = F (mod (5F — 4)) be soluble. 

If F = (mod 3), suppose (0, 6) gives an acceptable solution. Then working in 
Z3, from (|i:Ti)) = 2s + l, so s = 26'+l, from (irT^ s(s + l) = FA^-2A+l, so 
A={2e+ l){2e + 2) - 1 = 02 + 1, from (|4TT|) H ^ FA - 1 ^ 2; iio p = AH + 1 = 
29^ = 2 unless = 0, p 0; from gi!]) S = 2p - s = 6*2 - 261 - 1 = (6I - 1)^ + 1, so 
q = BH +l = 2{e -If = 2 unless 61 = 1, g = 0; with = 1, from JUl OOl [4T4|) 
C = B + 6l = 612-61 + 2, so r = Ci/ + l== 261(6' - 1) + 2 = 2 unless 6I = 2, r = 0. 
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Thus in Z3, for any 6 exactly one of p, q, r is zero, whence in Z 3\p, q or r and so n 
is only a C3N, possibly, if p = 3, in which case n — 561, with A = 1. 

For F = 5, (|4.15p is 0^ - 56*^ = -21, but = 5 (mod 21) has no solutions, so 
neither has (|4.15p . Thus there are no Ca^V's with A > G = rj = 1 and F < 7, and 
Theorem I4.3r b) follows in like manner to Theorem 14. 3f a) above. □ 

From (I4.14P for a solution to (|4.15p to yield a solution to (|4.13a|) we need (j> = —2 
(mod 2F) and 9 odd. Clearly from ([2T0| and ((49| with = 1, E ^ p-s ^ s'^+A > 
2, as required. Also such a solution will result in positive integers A, B, C, H, F 
which satisfy (|2.5ap . whence we have h :— gcd(A, B, C) = gcd(A, B) = gcd(A, C) = 
gcd(B,C); and from (g^ and s = 2 (mod h), whence from (|4.13ap 6=1 

(mod ft.), so /i = 1 or 5. Then if /i = 5, from (|4l4)) wc have (0, 9) = (3, 0) (mod 5); 
so if (0, 0) ^ (3, 0) (mod 5) then h = \ and our solution is acceptable. 

Using the theory of Pell's equation, ii (f>i'^ — F9i ^ = fc, a;^ — Fy^ = 1 and 

(4.16) = (2Fy2 + 1)0,, + 2Fa:y0, , - 2x#, + (2i^j/2 ^ ^^g^^ 

then ^i+i ^ — F9i^i^ — k and (/ii+i = (j>i (mod 2_F); so (|4.15p has an infinite 
sequence of acceptable solutions if it has a solution (^i, 6*1) with = — 2 (mod 2F) 
(unless, if possible, (^1, 6*1) = (3, 0) and y = (mod 5)). If F = s^ + s + 1, as for the 
ifa-family {712(5)} (equations (|3.4|) ) then (0i,6'i) = (2s(s + 1), 2s + 1) is acceptable 
and gives rise to 712(5) with A = 1 and B — 2p — S. Also there are acceptable 
solutions with F ^ s'^ + s + 1: for example for F = 87, {(j>i,9i) = (1912, 205) and 
yl= 11. 

For f 3, s = 1, so (01, 6*1) = (4,3), which gives 772(1) = 561; and for F = 7, 
s = 2, (01, 0i) - (12,5), gm is 0,+i = 1270, + 3360,, 9,+, = 480, + 1270,, 
and (02,02) = (3204, 1211) leads via gH ICTl Ol OOll to ^ = 229, s = 605, 
H = 1602, p = 366859, B = 733113, q = 1174447027, C = 734324, r = 1176387049 
and 2p-S - 733112.8118 <B = 733113 <2p-S+^- 733113.000737 for the 
K3N n — pqr. 

If 77(7) is the K^N arising from (0^, 9i) for F = 7, it is easily shown that, for i > 2, 
n{i + 1) — 254*77.(7) — 1.73 X 10^^77(7), and a naive "probability" estimate based 
on the knowledge that 77(2) is not a Cj,N and an assumption of the independence 
of the primality of p, q and r is : "p" (77(7) is a Cy,N for some i > 2)— 2^00- 

Challenges 2(a): Find a. C3N with B > 2p - \Jp~\ +1 + ^75- 

2(b): Find a C3.N with B > 2p - ^Jp~\+\■ 

From Theorem 13.21 r' = \p'{p + 1)^ is possible, so given p we can have 
C = \ [p + 1)^- Again, we show that this can be slightly exceeded: 

Theorem 4.4. For any K^N, C < ^{p"^ + 2p + ^^/4p~~S + with equality iff 
E = 2, F^G=1. 

Proof. For any K^N, A < p, and so from (|2.16p and (|2.19p we have 

But F > 2 and G > 1, so if F > 3 then C<4^+p+f = H^^^tll < 1 + 1)2 fo^. 
p > 3, and if G > 2 then G<^^ + f + f = SiE+^l < i(p + 1)2; thus G < i(p + 1)^ 
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unless E = 2,G = I in which case, smce G = 1, from (|4?T2l) FA^ - A - p' = 0, 



whence A = ^ + V4fp' + 1 ^ _L+ ^ ^4^. + i. < i(i+^4p-3), with equality 



iff F = 1. Hence from ((4T7| . C < ^ + p + ^{1 + V4p- 3), and Theorem 113 
follows. □ 

For a Ka-family with equality in Theorem [i^l from (P7T5|) . (P?Tg|) and we 

get 

(4.18) A := i7 + 1, B = AH + 3 = B{H) := + H + 3, 



C ^ C{H) :^ i + 2H^ + 5^^ + 

p = P(i/) + ij + L g = Q{H) := + H^ + 3H+1, 

1 
2 



5i/ 



:= i (h'-" + + 5i?-'' + bH^ + 4iJ + 2 ) and 



(4.19) 



2' 

With H — 2h this yields a iiTs-family with parameter /i, but from (|4.18p ii H = 2 
(mod 3) then gcd(A, B) = 3, so we take H = 6t ot H = 6i + 4 to get two 
iCa-families with the required maximal C property. Also if = 1 (mod 3) then 
3\p and 3|g, so for CsN's we must have the i^a-family with H — 6t. A search 
by Matthew Williams (see §5.ip up to i = 1365 found CsiV's only for t = 1 and 
t = 210, giving n = 43 • 271 • 5827 = 67902031, with ^ = 7, B = 45, C = 971, and 
n = 1588861 • 2001967381 • 1590423947471521 = 5058896665381789187674264635361, 
with A = 1261, B = 1588863, C = 1262241228152. 

The above _R"3-families will be shown to exemplify further bounds in Theorems 
UandSIOl 



4.2. An upper bound for q, given p and n. 

/Tl IP 
V 12' 

TL Tl 

Proof. Since r > q + we have q = — < — , whence pq^ + pHq — n < 

pr p[q + H) 



n H n H fTT 

and q < 1 / — \ . Now W — I < , / 

Vp4 2 \ p 4 2 y p 



4 y p 2 V 12 2 



iff J^ + ^-.I^<^-M^I(h- 



so Theorem 14.51 follows if we can show that — ^ < - ( H — \ ^ ) , i.e. 

SVbC 2V V 3 " 



20 



J. M. CHICK 



We consider two cases: (i) if -g- < then since S > 2, C > 3 and H >2, 




3iJ2 



> 1- 



4V6 



1 13 A 
12 18 SH' 

13H 



whence (|4.20p foUows. 



(ii) If f > ^, then n is a UgA-K^N, B > A > ^ and C = 3H + j > 3H (since 



7 > 0, see the proof of Theorem 14. li just after (|4.ip ). From Theorem 14.1 



A < m- 



H A 1 

so 1 

2 3H 3i?2 



<1 ^ + -^. Now l-(l-x)2>ix 



for < X < 1. Hence AVBCi l-\ 



1 



3H 3H^, 



>V(T)--{K^-3^)} 



1 \ 

---)>1 for i/>4, 



which is (gSQl); and for iJ = 2, since ^ < A<3H - y^-f, we have 4i < A < 5, 
which is impossible, estabhshing the result. □ 



For the C^N 191 • 421 • 431, Theorem |43] gives g = 421 < 421.981 



4.3. Upper bounds given n for p,A,B,C and ABC. A cursory glance at a 
list of C^N's suggests that a substantially better bound than the p < given by 
Theorem l3.4l should be attainable. That this is not so can be seen by considering the 
Chernick type ifa-families discussed in ^ 32. 21 with {A, B, C) — {2u — 1, 2m, 2u+ 1). 
For this family ^ 6u, = 12^^ -I, ABC = 2u[Au^ - 1), so ([23a|) requires 



H and F satisfying {12v?' 
F = 6u^- 5 and H = Ho = 
so the general solution is 



-1)H 



6u — 2u{Av? — 1)F. We see that for u > 1, 
1) is the unique solution with < Ho < ABC, 



(4.21) H^Ho 
In terms of i? = 



1 



Pf : 



If 



-tABC = 2u{2(u2 
2u, this gives 
Pf{u,t) 
qfiu,t) 



1) + (4^2 - l)t} 



2 ^ 



2(B^ - l)t}. 



= -B{B - 1){B^ - 4 + 2(5^ - l)t} + 1, 



= -B^{B^ - 4 + 2(^2 - l)t} + 1, 



= -B(B + 1){B' - 4 + 2(^2 _ i)t} + 1, 



and we have the two parameter system of i^s-families nf{u,t) :— PjQfrj for 
u — l,t > 1 and u > l,t > 0. If now we arbitrarily describe any KN with 
then nj:{u,t) gives flat K^N's for u > 3, since in general 

< 1.4 for B > 6; this also 



^ < 1.5 as "flat 

r ^ CH+l ^ C i i r ^ 



2 

B-1 



shows that for any Chernick type i^a-family the first member is the "flattest" , since 
"steepness" := ^ increases with p (and hence t). In (12J Pinch says that for CiV's 
up to 10^^ the largest value of p occurring is 72931, dividing 651693055693681 = 
72931 • 87517 • 102103; this is ri/(3,69). Gordon Davies did a cursory computer 
search which found the C^N^s n f{5, 3) , n f{8, 0) , n f{17 , 1) , n f{51, 0) and the very flat 
n/(102,0) = 861618073 • 865862497 • 870106921 = 649136982888522736355512801, 
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with "steepness"— 1.00985. Clearly any improvement on p < yfn will be of the 
form p < (1 — 0(1)) -^/n, with o(l) > 0, and we now show that 

Theorem 4.6. For any K^N , 



9 



Proof, li k := B — A and I := C — A, for any K3N we have n = p{p + kH) [p + IH) . 
From Theorem H?!! A < 2,H - I so 



(4.22) > > = -, 

^ ^ 3 3 3' 

and hence kH > H > and for / > 3, IH > 3y^, while if / = 2 we have 
n e {nf{u,t)}, and IH = 2H > 3^/^ is equivalent to 4ij2 > 3p = 3AH + 3, 
i.e. to H{4:H — 3{B — 1)} > 3, which is obvious from (|4.21[) : so in all cases 
IH > 3^1 y/3p. So for any K3N, n > p{p + ^/^){p + ^/Sp), and if n = 



Y + voV^), then p < V. Using the approach of Theorem 13. 51 to express 
V in terms of n as a power series, if we put n = 27y^, V — 3v^, v = y + J2^o '^iV~^ 
and equate coefficients of y^, y^ and y'^, we get Oo = — |, ai = j^t^2 = and 

^ = 32/^ - |y + 1 - rii + • ■ ■ . whence V = ^ - ^ ^ + I - + . . . But 

p < V and p is an integer, so Theorem 14.61 follows. □ 

From (|4.6p for n* we get p = 547 < \^ ~ = 574, and clearly a 

marginally lower bound than that given by Theorem 14.61 could readily be estab- 
lished. I speculate that for large n there is an upper bound y^ — fiy^{l + o{l)) for 
some /i > 0; such a bound can be established for the system nf{u,t) with fj, — 
But n* (see (|4.6p ) is both slightly flatter and much smaller than n/(5,0): if (using 
the notation of the proof of Theorem 14. ip there were an infinite system of bigA- 

K^N^s with (a,/3,7) = y "f-i 1: y -y); from (|2.6p we can deduce that it would 



require u = 4^— , suggesting a conjectural upper bound as in Challenge 3 below. 

V 108 

But possibly there lurk yet flatter K^N's capable of defeating any ^ > 0. 

/ n 5_ 6/^"_| 5_ 12 / n 

108 18 V 27 ~'~ 72 V 216 



Challenge 3: Find a K3N with p > {/^ - 4 4 '{/^\ , or prove 



that none exists. 

Theorem 4.7. For any K^N, 

A < VS^- 



2x02 



Proof. From Theorem 14. 2[ A < ^/3p — ^ \/^, and from Theorem 14.61 we have 



2 V 12 ' 

p < -y/n — + 1 < giving the result at once. □ 



Clearly a slightly tighter bound could easily be found. For n* (see 14. 6p this gives 
A = 39 < V3^ - - 40.81. 

If Bt(n) is to be a bound for aU K^N's for B given n for which inf (i?^(n) — _B) = 
0, then i?^(n) as in Theorem 14.81 (a) below suffices. But if we are content with 

inf (^) = 1 the n we could much more easily establish _B^(n) = \/2^/n. However, 
we outline a proof of 
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Theorem 4.8. (a) For any K3N, B < %/2^- (\/3 - ^ 
(h) For any C3N except 6601, B < - f V? 



39 + ^/3 
24 

1\ 175-3^/7 



2/ V 4 56 

Proof. We define B,in) ^2^- (v^- 1) ^+ ~ + and 

then (a) states B < B3{n) and (b) states B < Bj{n). Also for any KN, n > 561, 

/ TL 

and so Bjiji) < B^in). Further, if we put u — y we have 

(4.23) S,(n) ^ B;{u) := 2.^ - (^^7- i)u + ^ - ^ 

We note that for < 7, Bf^ and _B* are increasing functions for n > 561 (actually, 
for n < 62.7). 

We consider three cases: (i) G > 2, (ii) G = 1, F > 2 and (iii) G = = 1. 

(i) For G > 2, from Theorem [U and flJ^ . B < p. Also from > h 
so p'(ij2 + i?) > p'(£ + 1) = > 2^ for p > 3. Hence, with (EH) and since 
G > B, we have n > n' = ABCH{H^ + F) > B^p'{H^ + F) > > ^, so 
S < \/3n. This immediately gives B < \f2^fn, but the tighter bound B-iiii) then 
requires \/in < B-j{n), which holds for n > 5.625 x 10^; a computer check verified 
TheoremSHfor n < 5.625 x 10^. 

(ii) For G = 1 from (j4TT|) F is odd, so F > 3, and H ^ AF-l, so H"^ = Fp'-H. 
Also ?7 > 1 and from ((2?20)l ^ = ^ = H <p'. Hence, and from ((2TT5)) . we get 

n = p{H^{4p^ + rip) + H{bp + ?/ - 1) + 1} > p{i/^(4p2 +p) + ^Hp + 1} 

= p{{Fp' - H){4p^ +p) + 5pH + 1} = p{Fpp'{4:p + 1) - 4pp'i7 + 1} 

> p{3pp'{Ap + 1) - 4pp'^ + 1} = Sp"^ - p^ - 7p'^ + p 

So if we put w :— 'w{p) :— ^|(8p^ — p'^ — 7p^ + p) we have w < u — ^/^, 
and with ^ = 7 in (|4:23)) we get B^w) < B^{u) = B^{n) for p > 3. Also from 
Theorem l43l we have B < B{p) := 2p - \Jp^j + ^i^- We now want B{p) < 
Bj{w) for p> 3, i.e. /3(p) := Bj{w{p)) — B{p) > 0. The dominant term in P{p) is 
2(^/2— l)p, so /3(p) > for sufficiently large p, and /3(3) — 0.596 > 0; an outline of 
a general proof is as follows: using the substitution p = p{v) := v'^ + v + 1, we get 
from w{p) above 

w ~ (v) := w{p{v)) 

(4.24) 



32i;7 + 79i;6 + I25t;5 + 139t;4 + I07v3 + 54t;2 + 16w + 1 



and also bI{v) := B^{w^{v)), B^ := B^{v) := 2p{v)~^p{v) - | + i = 2i;2 + t; + 2, 
and then 

(3\v) -.^mv)) = B\{v) - B,„iv) - -i= 
=2u;2 _ B,n ~{V7- -)w + 



2' 56 2V3 



> — ^- > -Iv^ + 6v + —), 

5041)6 5V 4/ 

^t/ N 1/2 n 57\ / /- 1\ t/ N 4u2 + 24u + 57 / /- 1\ 
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Lemma 4A. If x > 0, y > and a;"* > y'^, then x > y and 

4 4 4 4 

a; — y a; — y 

^ x'^ + x^y + a;?/^ + Ax^ 

Applying this Lemma in (|4.25|) to 2w^ — Bm, with (|4.24p we get 

16t;8 + 96?;^ + 228^^ + 396w5 + 411?;'* + 324i;3 + 128i;2 + 32u - 12 
~4(32i;8 + 128i;7 + 316v^ + 500v^ + 556i;4 + 428^3 + 216i;2 + 64w + 4)i 
Now for w > 10 we have 50v^ > 16(u'^(w))*, so 

+ 24w^ + 57v^ 

and hence from (|4.25p 

57\ / r- 1\ . 4w2 + 24u + 57 1\ «/50 

— V 
16 

4^2 + 24w + 57 5 4u2 - 26t) + 57 

> ^t; 7:v ^ > for v > 10 

20 2 20 

For V = 10, p = 111, and a computer check verified Theorem 14. 81 for all K^N^s 

with 3<p< 109. 

(iii) Let Ug be a K^N with associated if-variables Ho,Ao,Po, etc, with Fq = 
Go = 1. Then from (|2l8l) = i?o + 1, so po ^ Hi + Ho + 1, from (|2^ 
Ao = AoGo = Ho + 1, from (|2.14ap Bo = 2po + and A,, < Po, so the 

conditions needed for the proof of Theorem 13.11 and case (ii) of Theorem 13.51 are 
met, and in a similar manner we define real variables n,p,A,B etc in terms of 
independent real variables H and rj. Then, keeping H = Ho, constant, we get from 
the above with (|2.13ap . (|2.14ap and (PTS)) 

(4.26) B = B{H, rj) := 2H^ + 2H + 2 + ^{t] - ^/At^H"^ + At]'{H + 1) + jf) 
and 

n = n{H, r/) :^ AH^+12H'^ + (24 + r^)H^ + (33 + 2r])H^ + (34 + 37])//"* 

^"^■^"^^ +(26 + 3??)ij3 + (14 + 2ry)ij2 + (5 + ri)H + 1 

Then as for Theorems 13 . 1 1 and 13. 51 we have that as rj increases, so B decreases and n 
increases, so if 1 < ?7a < ?7o, then B{Ho, rjo) < B{Ho, T]a) and n{Ho, rjo) > n{Ho, rja). 
So for the best possible bound we want to choose rja to be the smallest possible 
r]o (we already know from the proof of Theorem 14.31 that 770 = 1 = Fq = Go is 
impossible for K^N^s). Dropping the zero suffixes for our K^N, from (|2.12p we 
have s^ + T]s = ri{H'^ + H+l)-{H+l) = rjH'^ + -q'H + 77' and with (|2.13b|) we get 

(4.28) := 2-qH + 77 - 1 and 9 = 2s + rj and then 

(4.29) <f -riO"^ = -{if + iTf -2ri^l) 

Compare this with (|4.14p and (|4.15p and the accompanying discussion of Pellian 
solutions, which we now apply to (|4.28p and (|4.29p . We define acceptable solutions 
to (|4.29p in the same way as for (|4.15p . and for a solution to be admissible we require 
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9 — T] even, and </i = 77 — 1 (mod 4?/) since in ()4.28p we have H even. Henceforth 
replacing F with 77 in (I4.16P with — rjy^ ~ 1, if 77 is not a perfect square and 
(|4.29p has a solution, then (|4.16p gives an infinity of further solutions. We find 
that there are no solutions for r/ = 1, 2 or 5; for 77 = 4 there is the unique solution 
= 51,61 26 which leads to the KgiV 43-451-607 = 11771551 with B = 75; while 
for each 77 e {3, 6, 7} there are four fundamental solutions from which all other 
solutions can be derived via (|4.16p . but the only admissible ones are = 14, 6* = 9 
for 77 = 3 and (p = 90, 6' = 35 for 7; = 7, for which 0^ — 79^ = —475 reduces to 
X2 _ ^ _iQ via (l)^5X,9 = 5Y. 

With reference to (|4.16[) . amended as above, we need to show that for t] < 7 
any solution (cjji, 9i) is admissible iff its fundamental solution 9i) is admissible. 
From ()4.16|) 9i+i = 9i (mod 2), so (^?i+i — 77) is even iff {9i — rj) is even. Also the 
inverse transformation for ()4.16|) has (j>i — (277?/^ + l)<j>i+i — 2r]xy9i+i, and we have 
Vy^ — (2; ~ 1)(2; + 1), so if xy is odd then 8 | 7/, whence xy is even for rj < 7. 
Working now in with i] < 7, if 7/ is even then 4>i+i = 4>i, and if 77 is even then x 
is odd and y is even. But if odd rj = 21^ + I and y is odd, suppose (f)i = 77 — 1; then 
4>i+i = {2riy^ + 1)0; = (2?;+ 1)(77— 1) = 2iy{2i]+l) = Avrj + rj— 1 =77—1; conversely 
by the inverse transformation 0^ = {2riy^ + l)0i+i, so 4>i+i = 77 — 1 iff (/j^ = 77 — 1, 
whence by induction (jji = 77 — 1 iff 0i =77—1. Thus as required {(pi, 9i) is admissible 
iff (01, 9i) is admissible. Also for acceptibility we require E > 2: it is easily shown 
that E = F ^ G ^ 1 gives rise to 7; = ff^ + 2H^ + H'^ + H + 1 (cf (|4l8l) and 
Theorem 14. 4p . so 77 = 39 (with H ~ 2) is the least 77 with an admissible solution 
with E ^1; thus S > 2 for 77 < 7. 

With 77 = 3, the solution = 14, 6* = 9 gives the C3N 7 • 23 • 41 = 6601, but this 
is the only C^N, since we find in Z7 that the cycles given by (I4.16P are of period 4, 
and in Z with obvious notation gcd(^4i-|_2, i34i+2, (^41+2) = 7, not acceptable, and 
P4i+i = P4i+3 = qa = ra = (mod 7). This is sufficient for our proof, but for 
77 = 7 the solution = 90, 6* = 35 gives the C3N 43 • 433 • 643 = 11972017 with 
B — 72, and there seems to be the possibility of further CaiV's in the sequence of 
KsN's generated by dHH). 

It remains to show that for 77 = 3 or 7 and any even H > 2, B{H,ri) < 
B,^{n{H,Tj)). If we write (ji?^ as B = 2H^ + 2H + ^ - ^H{1 + x{H,7j))i and 
(|4?27l) as f = i?*(l + y{H, 77)), then for large H with x := x{H, 77) and y := y{H, ij) 
we have x = 0{jj), y = 0(i) and B^{n) = 2H\l + y)i - - + 7/)i + 

— Using the binomial series and omitting the complicated details we get 

i.,(7.)-i..f^-^ + iM + iiZ_^_^)l + ofi,)>o 

' V 32 64 128 256 4^ S-q^J H \H^J 

for sufficiently large H. To prove the result rigorously for all even H = Ho, we can 
truncate the various binomial series and approximate x and y to form functions 
bri(n{r], H)) = brj{n) and b{r],H) — b such that Bjj{n) > bri{n) and b > B, with 
briin) and b containing a relatively small number of terms all of which are retained, 
whence for each 77 we can determine a precise H*(rj) such that bri{n) — 6 > for 
H > H*{rf). My method of truncation and approximation was as arithmetically 
economical as I could make it, subject to retaining exactly the above term in jj, 
and after heavy detail arrived at H*{2>) — 240 and H*{7) ~ 66; Gordon Davies (see 
^5.ip did the computer verifications iov 2 < H < H*{7j). 
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Thus for any K^N with Fo = Go = 1, we have rjo > 3 and for any C^N except 
6601 we have t^o > 7, so with r^a = 3 or 7 as appropriate we have rjo > tja and 
Bo = B{Ho,Vo) < B{Ho,7]a) < B^^{n{Ho,Va)) < B^^{n{Ho,Vo)) = B^^rio) as 



We can construct a i^s-family {ri4(i)} such that, for any fixed /i > 1 and large t, 
B{t) - B^,{ni{t)) ~ B^(^i){ni{t)) - \/2-yn4(i) as foUows: = G = 1 and then 
H = t{t + \), A = -q = + t + I, s = tA, giving p4(t) = + 2t^ + 2t'^ + t + 1, 
B{t) = 2t4+3t3+3t2+t+2, C{t) = 2t4+5t3+6t2+4t+3, (?4(<) = 2t^ +bt^ +Qt'^ +A£^ + 
3^2 + 2t + 1 , r4 (t ) = + + 1 li"* + 1 Ot^ + Tt^ + 3i + 1 and ^4 ) = p4 (t) • 94 (i ) • r4 (i ) . 
Then we have B - 2f^ and n - 4ii6^ so B - B^(n) - B,,(n). Then 

n4(l) = 6601 and 714(2) — 11972017 as above, but we found no more GaiV's up to 
t = 31. Obviously (|4.29p is satisfied identically by the parametric forms for H,ri 
and s of {n4(i)}, since F — G — 1. 

For 7? = 7 and ((/'i,6'i) = (90,35), (OBll gives (02, 6*2) = (23190,8765) which 
yields the K^^N n = 2743993 • 9080853193 • 9095368033 with B = 5483607, H = 
1656 and B'j{n) = 5483607.001 (calculator accuracy). 

We note that for the sequence of K^N's associated with Theorem Upland (j4.16p 
for given F > 3, B ^ \/2^fy. I confidently conjecture from the above proof of 
Theorem 14.81 and from numerical evidence that B < Bri{n) for every K3N, but 
have not attempted a general proof. 

Challenge 4: Find a C3N for which n > 6601 and Bjin) - B < 0.1. 

To conclude 2] we shall show that the K3- families deriving from ()4.18p and 
()4.19p . with E = 2 and F = G = 1, not only give equality for the upper bound for 
G given p, as described in the discussion following Theorem 14.41 but in the same 
way give equality for upper bounds for G, BC, ABC and n given H , and for G and 
ABC given n. Using the notation of (|4.18p and (|4.19p : 

Theorem 4.9. For any K3N , (a) C < C{H), (b) BC < B{H)-C{H), (c) ABC < 
A{H) ■ B{H) ■ C{H) and (d) n < N(H), with equality iff E ^ 2 and F = G = I. If 
there is equality then H ^ 2 (mod 6), and if also n is a C3N , then H = (mod 6). 

Proof From ([2T8)) . A = ^i^, and then via (|2.19p . (|2.16p and (j2.6p we can express 
B,C,n ( and obviously also p,q,r) in terms of E, F,C, H. We now regard all 
other variables as functions of continuous real variables E, F, G, H, subject to all 
the relations so far established for K^N's. Let Z be any of, or any product from, 
A,B,C and n': regarding H as fixed, we write (j)z{E, F,G) :~ Z evaluated at 
{E, F, G). From (HI]) and (i2T8)l we have n' = ABCH^ + BCH{G + H), and then 
we readily see that F occurs only in the denominator of any Z, and thence 

(4.30) (j)ziE, F, G) < (j)z{E, 1, G) unless F = 1 



Lemma 4B. // f{x) := ax ^ — with a > 0,b > 0, then for < Xi < \ — we have 



required. 



□ 



b 



X 




/(xi) > f{x) iffxi <x< . 



This follows at once from 



f{xi) - f{x) = a{xi - x) + b( ) = (x - xi)( a 



\Xi X/ \XXi 
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From ((2A61 [2l8l [219]) wc have 

ABH + A + B 1 uH + G\fH^ + GH E+1 



F J\ FG G 
H + G /H^ + GH E + 



H 



FG G 



I.e. 



(4.31b) G ^ <t>c{E^ F, G) ^ ^{(g^ + F)G + } 



2ij3 H 3H 1 1 



£;^2 F EF G EG FG 

Obviously since E > 2, 

(4.32) (/)c(£;,-F,G) < (/)c(2,F,G) unless £; = 2 

Similarly, the last three terms of (j4.31bp increase as G decreases. Also, applying 
Lemma 4B with x = G, f(G) [H"^ + F)G + and xi = 1 since G > 1, 

we have /(I) > /(G) iff 1 < G < ^'(^'^+2^) . But from Theorem [U G < 2H, 
and since H > 2 we have 1 < G < 2H < < ^^( ^l+^J ), so /(I) > /(G) unless 
G = 1. Hence from (|4.31bp . 

(4.33) (/)c(£;, i^, G) < (jjciE, F, 1) unless G = 1 
Applying (|430l 1021 14:33)) in turn, 

(4.34) G = 0c(^^,^,G) < 0c(2, 1,1) unless £; = 2 and F = G = 1 
Combining this with the discussion of (|4.18p and (|4.19p . we have Theorem 14. 9f aV 

Also B = § + ^ + ^ (as in (lOTal) ). so (j)B{E,F,G) < (j)BiE,F,l) unless 
G = 1, and hence with 

(4.35) BC ^(j}Bc{E,F,G) <<j)Bc{E,F,l) unless G=l 

Now given F^ G and iJ, also determined are A — and p — AH + 1. Hence 

from ((2T91) and ((2T6)) we have 

+ E\ fBp + A\ p + E fp'^ +pE + AG\ 



^Bc{E, F, G) = j [-^) =^[- 

G2l 



EG J 

^ P^ + AG^ 2p^ + AG 



G2l E J G2 

Applying Lemma 4B with x ^ E, f{E) = E + and xi = 2 since E > 2,we 

have /(2) > f{E) iS 2 < E < which is true unless E = 2, since p > 3 and 

by Theorem [121 E <p-l. Thus 

(4.36) ^BciE,F,G) <(j)Bci2,F,G) unless = 2 
Hence from KM 11351 , 

(4.37) BC ^ (t)Bc{E,F,G) < (l3Bc{2,lA) unless = 2 and = G = 1, 

and Theorem I4.9r b) follows as for H.Qi a.) . 
Also A is independent of so from (|4.36p 

(4.38) (pABciE, F, G) < (l)ABc{2, F, G) unless E = 2 
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R,rther, AB = p. = " + + ''^ + '''^ + '> 



F \ FG 

+ F{E + 1) \ 2H^ + F{E + 1) 



G I F2 



so applying Lemma 4B with x = G, /(G) = G + ^ +-^^-^+1) g^j-^^j x\ — \ since 
G > 1, we have /(I) > /(G) iff 1 < G < i/^ + + l), which is true unless 
G = 1, since i? > 2 and from Theorem [131 G <2H. Hence 

(4.39) (j)AB{E,F,G) < (I)ab{E,F,1) unless G = 1. 
Then from (|4:33l) and gJH), 

(4.40) (l)ABc{E, F, G) < cj)ABciE, F, 1) unless G = 1, 
and from KM KM ODD 

(4.41) ABC = (jjABciE, F, G) < ?!)asc(2, 1, 1) unless = 2 and F = G = 1 

and Theorem 14. Of c) follows as for l4.9f aV 

Also from and g3Dl), (j)ABc{E, F, G) < 0yisc(2, F, 1) unless 2, G = 1 

and from (EH) n' ^ ABGH{H^ + F), so 

^n' {E, F, G) = 4>ABc{E, F, G)H{H^ + F) < (j)ABci2, F, l)H{H^ + F) - (2, F, 1) 

unless F 2, G = 1; and from ^^M) 0„'(2, F, 1) < 0„-(2, 1, 1) unless F = 1. Thus 
n' = 4)n,{E,F,G) < 0„'(2,1,1) unless F = 2 and F = G = 1; Theorem gH follows 
as forllija). □ 

We now express our results for upper bounds for G and ABC given n in terms 
of the functions of (|4.18p and the inverse function of (|4.19p : using the method 
of Theorem 13.51 we could also express our results as series of descending powers of 
v^, but we simply indicate the leading terms: 

Theorem 4.10. For any K3N, (a) C < G(iV"i(n)) and 

(b) ABC < A{N-^{n)) ■ B{N-^{n)) ■ C{N-^{n)) with equality as m Theorem\J^ 
For large n, C < 2~ini (1 + o(l)), ABC < 2-T^n^T(l + o(l)). 



Proof. Suppose that all variables are real as in the proof of Theorem l4.91 and that 
{E,F,G) = (2,1,1); then n = N{H),H = N~^{n), so A = A{N~^{n)), etc. Also 
the first part of the proof of Theorem 13 . 6 1 applies . with i = d = 3 and Xd = E = 2, 
giving r = y-Mll±l+l_ B^j^-t G = ^, so we have the functional relationship 



(4.42) C{N-\n)) = ^ 



4iV-i(n) 
Similarly from (|2.6p we get 

(4.43) A{N~^n)) ■ B{N-\n)) ■ C{N-^{n)) 



iV-i(n){(iV-i(n))2 + l} 

Now suppose that n is any K3N with its standard A to H integer set. Then 
(a) G = ^, from Theorem O r' < ^^"+^~^ with equality iff F = 2 and from 
TheoremllHi/ > N-^{n) with equality iff (F, F, G) = (2, 1, 1), and hence 
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C < "^^(^^ = C{N-'^{n)) from filiS)) . with equahty iff {E, F, G) = (2, 1, 1). Also 
(b) from (121]) 



ABC 



7? ' 77 

from (HIS]), with equahty iff {E,F,G) = (2,1,1). The first sentence of Theo- 
rem 14.101 foUows. 

Also for large n and equality, from (I4.18P and (|4.19p we have C - ^H*, ABC - iff^ 
and n ~ so H ^ {2n)T« and C ^ 2~in^ , ABC ^ 2~'^ni' , giving the final 

part of Theorem 14.101 □ 

5. The algorithm and its implementation 

5.1. Acknowledgments. This is the most convenient place to acknowledge my 
immense debt to my two friends who have carried out the computer implementa- 
tion of my algorithms: firstly Gordon Davies, like me a retired teacher of math- 
ematics at Haileybury College, England, who did the computing throughout the 
development stage, programming in BASIC V with 32-bit arithmetic, and by Au- 
gust 1999 taking us on his RISC-PC, running RISC OS 3.7 with 16 Mb of RAM, to 
C3 (2 X 10^*) = 42720 (where CaiX) is the number of three-prime Carmichael num- 
bers up to X); and secondly Matthew Williams, a recent Haileybury student and 
Cambridge University computer science graduate, who then got us up to C3(10'^'*), 
using his 1 GHz Athlon with 900 Mb of RAM and programming in C++ mostly 
with 64-bit arithmetic. 



5.2. Notation. In this section we shall use the following upper bounds based on 
results in § [3] and § [4j Each of these could be replaced with a slightly greater and 
simpler bound with only marginal loss of efficiency. We seek all CsN's less than 
X, where for convenience X is not a C3N, so n < X . 



fx-'-^^ 

9 



(a) p < pm := 

(b) A < Am 



ic)B<BM:=2p-^jp-l + ^^ 




(d) <7<Qi I 2P 

(e) q<Q2- 




(Theorem 

(Theorem 113) 

(Theorem 

1 (Theorem [3+ 



(Theorem 113 

p V 12 
(f) Z:=min(Qi,Q2 



(g) For any g < 
rem 13^ 




(Theo- 



5.3. A brief description of four algorithms and the split range procedure. 
In [12] Pinch describes the two algorithms which he used to find all CiV's up to 
10^^. We briefly describe these now, since we used them, slightly modified to take 
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advantage of d = 3, to do selective checks on our results for large X. In all these 
algorithms, for given X the outermost loop runs through all odd primes up to pM- 
In Pinch's first algorithm, which as modified by me for c? = 3 we call PI, for each 
p, E runs through the range El{p, X) < E < p' f Theorem I2.2p . where El{p, X) is 
a fairly complicated function, not given here, which I formulated using d = 3 (so 
P = p < q) to cut out some of the cases which would result inn > X (for p < \/3X, 
El{p , X) — 2); using El{p, X) reduced the time for PI by about a third. For each 
E a range of integer values of D is found, subject to 1 < A = DE — p^ < 2p' 
(Theorem 12.31) . and for each {E,D) pair q and r are calculated from (|2.9p . E, and 

D for each E descend through their ranges, and if r > — next E is taken; else q, r 

pq 

and Ai = (iHH) must be integers, with next D at the failure of any test, 

P' 

and q and r are tested for primality. 

For large X the remaining algorithms are all significantly speeded up by the 
split-range procedure, which we briefly describe. Suppose that variables x and y 
are connected by the bilinear relation axy + bx + cy + d = with a > 0, V := 

be — ad > 0, and that xi > , so over the interval xi < x < X2, y decreases 

a 

as X increases; and also that we wish to find integer pairs {x, y) over this interval, 
and that a trial where we start with x e Z (an rc-trial) costs k times the cost of 
dy 

a w-trial. Then if — = — fc at (£,,'4') we minimise the cost by using a;-trials for 
dx 

X < and y-trials for y < (so x > ^). So if a;i < ^ < a;2 it pays to split the 

range at [£,, tp) ~ I ^ , (as may easily be shown), with a saving, 

V Vka a ) 

compared with using only x-trials and ignoring the cost of deciding whether to split, 

^ . , ka(x2 — . , 

of approximately y-trial costs. 

axi + c 

Pinch's second algorithm as modified by me for d = 3 (PII) for each p runs 
through all primes q satisfying p < q < Z . For each (p, q) pair it uses the Euclidean 

algorithm to find H and hence A ^ and ii A > Am it takes next q. Else 
p'q' 

it finds Li :— = lcm[p', g'], i? = pq, and w such that wR = 1 (mod Li), 

H 

by the reverse Euclidean algorithm; then, since n — rR = 1 (mod Li), we have 
r = w + uLi, and also R' = Er' (|2.1[) . so eliminating r we seek integer pairs (u, E) 
such that 

(5.1) LiEu + w'E - R' ^0. 

With u ascending and E descending we use the split range procedure, take next q 

when r > — , and for each integer pair (u, E) we test r for primality. 

Our first successful algorithm (HI, originally devised when seeking Perrin pseudo- 
primes, before we knew of other algorithms) is the same as PII as far as finding 
A < Am- It then found r in essentially the same way as our main algorithm HII, 
described next. 

HII was motivated by the realisation that as p becomes larger in HI many more 
pairs {p, q) will result in H small enough to give A > Am ; and that by first analysing 

p' 

p' = AH such pairs need never be considered. Since H is even, A divides — and, 
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taking the most favourable case as an example, if — is prime (and so is a Sophie 

Germain prime) then A = 1, H = p' is the only possibility and all possible values 
of q will belong to the arithmetic progression (AP) q = p (mod p'), with q < Z. 
So for each p we factorise ^ to find all possible pairs {A, H) with A < Am, and we 
organise a set of AP's which will contain without repetition the resulting g-values, 

q' 

which we then test for primality. For each prime q we find B — we find integer 

H 

pairs {C\F) from (|2.5ap . which is bilinear in C and F, determining the range of F 
as described below in § I5.4i and test r — CH + 1 for primality. 



5.4. The implementation of the HII algorithm in more detail. The basic 
idea of the HII algorithm is as stated above, and we now describe more fully a 
method of programming it, drawing attention to certain worthwhile economies (our 
own method was slightly more complicated, as we shall explain briefly in the next 
subsection). We structure our description in terms of several loops, beginning with 
the outermost as Loop 1. To find all C^N's with n < X, we begin by preparing a 
bitmap prime database up to at least Q3 and we calculate pm- 

Loop 1: for each p satisfying 3 < p < pm, we calculate AmtBm and Q2, and 

p' 

form an array {yl(j)} of possible A values. To do this we have A \ — , so suppose 

P' 

that the prime factorisation = nfc=i Pfe''*"; where % > 1. We have A{\) — 1, 

and for j > 2 we can obtain A{j) by a process involving successive trial division of 
p' . . 

— by ascending primes, for each new pk multiplying all A{i) already found by pk 

to form more possible A values (but taking care to avoid duplication when 77^ > 2), 
testing for A < Am before adjoining A to the array, and storing the successive 
prime factors {p} {pi,P2, ■ ■ ■ , ps} of {A{j)} as they arise (so 5 < p). 

Loop 2: for each A{j) = A, we now develop a set of AP's which must contain 
q for any C3N associated with (p. A), with each AP having common difference p', 
and we also obtain the corresponding set of AP's with common difference A which 
contains the associated B values. We denote the A*'' term of the z*'' AP for q by 
q\{i), so qx{i) = qi{i) + A'p', and similarly for Bx[i). We have H — and for 
use in Loop 3 we define Fo = [f J and f := {Fo + 1)A - H. Clearly S = ^ < % 
and B < Bm; also in loop 3 we shall show that B < so an upper bound for B 
is (3 := min(i?7\/, %-,^)- Next we use the facts that q is prime and B is coprime 
to A to restrict the number of cases to be considered. For A > 1 we use a sieving 
method with those pk S {p} which divide A to find ^{A) :— {t : 1 < t < A' and 
gcd[^, t] = 1} and we define $(1) := {1}. Then for each t we form Bi = A + t and 
qi = BiH+l = AH+tH+1, and for ^ > 3, if 7? ^ (mod pk)^t€ ^(A) such that 
t = —jj (mod Pk) and then qi = Q (mod p^), so q\ = qi + A'p' = 9i + ^AH = 
(mod Pk) and thus the AP {q\} is entirely composite. Therefore if pk \ qi we do 
not adjoin Bi or qi to the arrays {Bi{i)} or {qi{i)} (Gordon found that up to 
X — 10^* this eliminated from consideration about 19% of the potential AP's). 
So we form the arrays {Bi{i)} and {qi{i)} of first terms and we use the iterations 
Bx{i) = Bx-i{i) +A and q\{i) = q\-i{i) +p' to form arrays {Bx{i)} and {qxit)}, 
having first examined each Bx-i{i) for associated C^N^s as described in Loop 3 
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below. Having used a sieving method to find ^{A), we get B\{i) and qx{i) increasing 
steadily throughout this process, and B\{i) > (} triggers next A{j). 

Loop 3: We write B\^i{i) = B and q\-i{i) = q, and if q is composite (check 
against bitmap prime database), we take next B and q. Writing K :— AB, U := 
KH + A + B = Bp + A and V -.^^ {A + B)H + 1 = p + g ', ((23a|) may be written 

(5.2) KCF ~VC-U ^0 

and we use this bilinear relation in C and F to find integer pairs (C, F) and hence 
possible r = CH + 1. Theorems 12.11 and 14.41 suggest that we consider the F-range, 
and we examine certain related economies including a procedure for splitting the 
range. 

Put r := ^ - 1 and T B + 1. Then r = ^ < ^ , so r' ^ C H < Y , and 
pq pg pq' ' 

usingdEll) we get F > /l := "^ky^ \ and C > T whence F < ^^W^- if 

Fl '■— T/lI and Fm '■— L/mJ , we need integer pairs (C, F) such that F^ < F < Fm- 

We easily show that 

, , P , pH{Bp + A){HB + l) 

fUB) -.-h -^ + J^+ ABiX-p-pHB) 

and /m(B) :=/m ^ ^ + + ^) + ^^(^ 

Clearly /l < Jm, but we can get Fm = L/mJ < /l < Jm < IIl] = Fl, in which 
case there are no possible F values. Further, as i? ti both fL{B) [ and Jm^B) J,, 
and also Fo = [-jj < < /l < ^l, so if Fm — Fo, then we take next A{j) (at 
X = 10^* Matthew found that this Fm — Fo trigger reduced the program time 
by 20%). It is easily shown that 

whence Fm(^) := L/m(^)J = -Fo, justifying B < anticipated in Loop 2. 

A further though smaller economy can be achieved by eliminating from consider- 
ation some or all of those B values for which Fo + 1 < /l < Jm < Fo + 2 and so 
Fl > Fm- if then fL{P) > Fo + 1, we can take next A; if not, it can be shown that 

if a := - *rqi£ ^^len fL{-^) - Fo + l, so we can jump to B - ^ to find 
/l(-B) just greater than Fo + l and then continue (but this is awkward to program). 

We next consider splitting the range (see and a method of economising on 
C-trials which arranges them in an AP, first term C'o, say, and common difference 
e G {1,2,3,6}. Consider the conditions (a) 2 | AB and (b) 3 | AS and 3 f H. If 
only (a) holds, C is odd and e = 2. If only (b) holds, for C = — -g- (mod 3) we 
have r = CH + 1 = (mod 3), so C ^ (mod 3) and C ^ (mod 3), leaving 
only one possible residue, and e = 3. If both (a) and (b) hold, then e = 6, and 
if neither, e = 1. Then by eliminating as appropriate for each situation over the 
range B + 1 < C < B + e if 2 | C, 3 | C or 3 | r, we find C. With the notation of 

and with a; = C, it seems reasonable to take k = - and then ^ = \/ ^ — VeH. 



K 

So to execute the loop, as described above, we see whether Fm and Fl values 
permit us to take next A (or possibly to jump some i3's); and then if Fl > Fm 
we take next B. li B < \J eH we do C-trials until C > V eH, taking next 
B li F < Fl while C < VeH, and then F-trials; but if i? > \/eH we do 
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F-trials for Fl < F < Fm- For C-trials F = where Ft U + VC and 
Fb :— KC, so with V* = eV and K* ~ eK we start with C — Co and then 
do Ft — > Ft + V* and Fb — > Fb + K* to find F for next C; and similarly for 
F-trials with C = j^p_Y — ^ we do E — > E + K for unit increase in F. Also if 
Fl gives _E = 1, we take next E. 

Each (C, F) integer pair then gives r = CH + 1 which we test for primality, using 
the standard algorithm if r is beyond the bitmap prime data base. 

5.5. Some notes on our implementation of HII. (i) In the development stage, 
to test q for primality we used a carefully designed but complicated system of 
tracking through a prime database, exploiting the advance of the arrays{gA(j)} by 
p' for each unit increment in A, and we also had much less RAM. For these reasons 
we constructed arrays {q\{i)} for each p, rather than each A, which involved extra 
complications with certain loop exits. But then Matthew found that primality 
testing for q was taking at least 80% of the time, and constructed the bitmap 
database, which ai X = 10^*, for example, reduced the program running time by a 
factor of at least 5, and was a major contribution to what we were able to achieve. 
Nevertheless we did not revise the array structure, as we estimated this would have 
given only a marginal decrease in time. 

(ii) In ^5.71 we shall give some running times , so we mention that the Loop 3 
Fm = Fo trigger was a late discovery, and right up to C3(10^'*) our implementation 
only used the special case Fm = 0, which Matthew's later trial showed gives about 
60% of the 20% time saving available at X = 10^^. 

5.6. A faster algorithm? If (C*, F*) is an integer pair, it follows from the theory 
of PII outlined in 15.31 (or directly from (|5.2|) ') that a necessary condition for (C, F) 
to be an integer pair is C = C* + Ku, and then (j5.ip and (|5.2p give 

(5.3) KuF + C*F~-Vu~C*F* ^0, 

an even more discriminating bilinear relation, between F and u. The total number 

of trials when the split range is used with ()5.2p is approximately 2\/H — T — —, so 

when this is very large (big H, very big X), using (|5.3p might be worth the cost of 
finding (C*,F*) — either by the method of PII for w, or simply using (|5.2p until 
(and if) such a {C*,F*) is encountered. We did not implement this. 

5.7. Comparison of algorithms for d — 3. Ignoring time required for primality 
testing of q in PI and HII (by virtue of "bitmap") and of r (relatively seldom 
required and the same for all four programs), and based on the number of test pairs 
{E,D),{p,q),{C,F) involved, I deduced that PI, PII and HI are all 0(X§+°(i)); 
in the range 10^^ < X < 10^^ for all three programs when X was multiplied by 10 
the multiplier for the time was close to 4.325 and slowly increasing with X, giving 
some support to this deduction since lOa — 4.64. On the same basis I conjecture 
but have not succeeded in establishing that HII is 0(Ar2+°(i)), ^^^^ ^^^^^ ^j^g range 
10""^^ < X < 10^* the corresponding time multiplier was 3.013, slightly less than 
105 ^ 3.162. 

Here are a few of the many times recorded. Illustrating the effects of improve- 
ments in computer technology, more powerful algorithms and the use of a compiled 
in place of an interpreted language, Gordon first reached C3(10^^) = 1000 early in 
1998 with an early version of HI on a mid-1980's computer in about 45 hours; in 
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November 2001 it took Matthew 0.19 seconds actual calculating "user" time (1.15 
seconds total) . Gordon on his new computer (see §5.1|) , with HII and my q primal- 
ity testing method (see subsection 15.5( 1)) in mid 1999 took about 32 seconds for 
C3(10i2) and 35i hours for C^ilO^^); in mid 2000 CsilO^^) took Matthew 10| min- 
utes with a slightly improved prime testing method and compiler optimisation, and 
finally with this method fully replaced by bitmap it took just 2 minutes 7.59 sec- 
onds. In July 2002 C3(10^"') took about 58 hours, with about 9 minutes for the 
bitmap. 

RAM and time constraints prevented us from going on to X = 10^*^. 

5.8. Checking and correction. Up to X — 10^*, Gordon and I had Richard 
Pinch's paper and his Internet results to check against. We achieved agree- 
ment up to lOi'^, but at 10^8 we found that his hst omitted = 835327 • 893359 • 
1117117 = 833645090806507981 (for more on nt, see discussion following Theo- 
rem 4.1). Richard told me that inexplicably failed to reach the Internet list, 
although his program gave it. He also kindly put me in touch with Carl Pomer- 
ance, who sent me the first preprint of [TU] and invited Gordon Davies and me to 
attempt the awkward evaluation of the constant K3 (see [B]). Some months later 
when Carl asked us for any counts we had beyond X = 10^^, Matthew had got to 
X = 10^°, but had not yet done any checks; it later emerged that a problem in the 
program was hy X = 10^° unfortunately causing omissions: the value of 120459 for 
C3(10^*') which we gave to Carl and is pubhshed in [10] should be 120625, and the 
number of imprimitive C3iV's up to 10^° should be 89854. 

Obviously comprehensive checking of HII results for large X with other known 
algorithms is not practicable. Soon after successfully programming PI, Matthew 
used it for a complete check a.t X = 10^^; this took about 62^ hours, checking 
q for primality by the standard algorithm, and no discrepancy was found. For 
final checking he used PII to find the C^N^s corresponding to every fc*'* p-value 
for X = 10^ for {N,k) =(19,2), (20,10), (21,30), (22,150), (23,1000) and (24,1000), 
with initial p values chosen to cut down repetitions of the same check. No discrep- 
ancies were found. The last of these checking runs, at X = lO^-^, took PII about 
74 hours and HII about 3-^ minutes. 

We are grateful to Harvey Dubner for collaboration which gave a further partial 
check. Let Cl{X) := #{n : n is a C3N with A = 1, and n < X}. In [8] Dubner 

finds C;J(10^) up to TV = 20, and suggests that for a "wide range of N" , J^) - 

63(10 ) 

0.644. He uses an entirely different algorithm for Cl{X), based on relevant (1, B, C) 
values. In correspondence he then took 6*3(10^) up to = 24, obtaining agreement 
with counts we have extracted from our discs for 6*3(10^^) and, later, 6*3(10^''). In 
Table 1 of g6] we extend up to = 24 Dubner's Table 2 for (1, B, C) in 0. 

When finding C3(A) for X > 10^^, we avoided the danger of rounding errors 
wrongly including or excluding a C3A^ very close to X by doing a run to find C3(A*) 
with X* = (1 -t- e)X and examining individually any C^N^s in the range (1 ± e)X, 
where typically e = 10"^ or 10""^ (at X = 10^" Matthew took e = 0.1). 

6. Statistics 

In Table 1 we tabulate for X = 10^, with 3 < iV < 24, C3(A) and various other 
numbers which we now define. In [lOj Granville and Pomerance define primitive 
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CiV's, and for CaN's their definition implies that a C3N is primitive iff H < ABC; 
C^{X) := =ff{n : n is a primitive C3-/V and n < X}, and our data are consistent 

C* (X) 

with their conjecture that — — - as X 00. 

CsiX) 

Let C := {n : n — pqr is a C3N and p = q = r = —1 (mod 4)}; Rabin showed in 
[15j that the probabihty of any odd composite n passing the strong pseudoprime 
test for a randomly chosen base b is less than j , and that this bound is approached 
most closely when 71 e C; and Pinch lists various other properties of C in [T^ : 
C{X) #{n : neC and n<X} 

of UQ] Granville and Pomerance conjecture that CsiX) ^ n— — ^ 

' J a\ J ^(logX)3 

, where T3 — 2100 is a constant whose evaluation is discussed in [6]; 

dt 



In 



27 



dt 



{\ogtr 

they define /3 and 7 by C3{X) 



X^ 



1_ 
27 



X3 



and predict that j3 



{logXf 27 J 2 (logi)3 
and 7 eventually converge to from above and below respectively. Our new data 
are consistent with this, supporting their cautious comment in [10] (but see [5], 
Table 3 and comment). 



Cl{X) is defined above in § 15.8 



Table 1. 



N 


CsiX) 


cux) 


C^{X) 
CsiX) 


CiX) 


/3 


7 


c'Ux) 


cm 

CsiX) 


3 


1 


1 


1 





32.96 


9.092 


1 


1 


4 


7 


7 


1 


1 


253.9 


53.13 


6 


0.8571 


5 


12 


12 


1 


1 


394.5 


78.07 


11 


0.9167 


6 


23 


19 


0.826 


1 


606.5 


128.1 


18 


0.7826 


7 


47 


36 


0.766 


4 


913.5 


220.2 


36 


0.7660 


8 


84 


59 


0.702 


8 


1131 


321.9 


59 


0.7024 


9 


172 


113 


0.657 


15 


1531 


519.9 


122 


0.7093 


10 


335 


208 


0.621 


29 


1898 


761.3 


227 


0.6776 


11 


590 


338 


0.573 


50 


2065 


961.5 


403 


0.6831 


12 


1000 


529 


0.529 


79 


2110 


1113 


680 


0.68 


13 


1858 


930 


0.501 


153 


2313 


1349 


1220 


0.6566 


14 


3284 


1550 


0.472 


271 


2370 


1496 


2104 


0.6407 


15 


6083 


2621 


0.431 


487 


2506 


1680 


3911 


0.6429 


16 


10816 


4201 


0.388 


868 


2510 


1763 


6948 


0.6424 


17 


19539 


6814 


0.349 


1569 


2525 


1839 


12599 


0.6448 


18 


35586 


11190 


0.314 


2837 


2534 


1899 


22920 


0.6441 


19 


65309 


18432 


0.282 


5158 


2538 


1947 


41997 


0.6431 


20 


120625 


30771 


0.255 


9443 


2538 


1984 


77413 


0.6418 


21 


224763 


51432 


0.229 


17316 


2541 


2019 


144300 


0.6420 


22 


420658 


85921 


0.204 


32351 


2538 


2047 


270295 


0.6426 


23 


790885 


143620 


0.182 


61130 


2531 


2067 


508780 


0.6433 


24 


1494738 


241562 


0.162 


115606 


2523 


2081 


961392 


0.6432 
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From Theorem ESI we have n < NsiP) := \{P^ + 2P^ ~ - + 2P^ - P). 

For any odd prime p we define x{p) ■= • nis a, C3N and n — pgr} and 

T{x) := J2p<xX{p)- Our hst of CsN's up to 10^^ enables us to find x{p) and 
T{p) up to p~= 11213, since we have A^3(11213) < lO^^ < 7V3(ii239). In Table 2 
we tabulate p,x{p) ^J^d T{p) up to p = 211, and in Table 3 p,T{p) for n{p) in 
intervals of 50 or 240 up to 7r(11213) = 1357. 

X(p ) = for p = 11, 197, 1223, 1487, 4007, 4547, 7823, 9833, 9839 and 10259, and 
x{p) = 1 for 51 values of p. x(211) = 17 is the greatest value of x{p) until 
p = 1171; with x{p) > 22 we have x{p) = 22 for p = 1171,7481,8521,8647 and 
10711, x(9241) = 24,x(10837) = 25 and x(2221) = 29. 

At first in Tables 2 and 3 T{p ) keeps remarkably close to p before going ahead 
for a bit, but then p gradually overhauls T{p) and seems to be slowly pulling 
away. Clearly C3{X) < T{pm), and a plausible heuristic argument that T{pm) < 
can be based on the loops of algorithm HIT, ignoring the primality 
requirement on p,q,r. The best upper bound for C3{X) which has so far been 
proved is 0(Xt4+°(^)), by Balasubramanian and Nagaraj in [T]. 

Table 2. Number and cumulative total of C^N's with first prime p 



p 


3 


5 


7 


11 


13 


17 


19 


23 


29 


31 


37 


41 


x{p) 


1 


3 


6 





5 


2 


2 


1 


2 


7 


5 


7 


Tip) 


1 


4 


10 


10 


15 


17 


19 


20 


22 


29 


34 


41 


P 


43 


47 


53 


59 


61 


67 


71 


73 


79 


83 


89 


97 


x{p) 


11 


3 


3 


1 


10 


3 


7 


4 


1 


2 


5 


6 


Tip) 


52 


55 


58 


59 


69 


72 


79 


83 


84 


86 


91 


97 


P 


101 


103 


107 


109 


113 


127 


131 


137 


139 


149 


151 


157 


x{p) 


2 


5 


3 


10 


5 


5 


11 


4 


6 


2 


9 


11 


Tip) 


99 


104 


107 


117 


122 


127 


138 


142 


148 


150 


159 


170 


P 


163 


167 


173 


179 


181 


191 


193 


197 


199 


211 






x{p) 


7 


2 


3 


4 


11 


6 


10 





7 


17 






Tip) 


177 


179 


182 


186 


197 


203 


213 


213 


220 


237 







Table 3. Cumulative total in Table 2 extended 



n{p) 


47 


97 


147 


197 


247 


297 


347 


397 


637 


877 


1117 


1357 


P 


211 


509 


853 


1201 


1567 


1951 


2341 


2719 


4723 


6823 


8999 


11213 


Tip) 


237 


565 


896 


1235 


1556 


1906 


2299 


2651 


4347 


6110 


7945 


9608 



From (Prg)) and Theorem HUd) we have n < N{H) = \{H^^ + 4iJ^ + lAH^ + 
30i/^ + bZH^ + 69i?5 + IIH'^ + 55i/3 + aii?^ + 12H + 2). 

Since A^(268) < lO^** < 7V(270), we can in a similar way use our list of C^N's up 
to lO^"* to count all the CaiV's for each iJ up to = 268. Let 

C{H) := #{n : n is a CgiV with gcd[p',<7'] = H} and Z{x) Eh<^CW. 
We find C{H) = for if = 68, 76, 160, 176, 188, 196 and 218, ((H) = 1 
for H = 98, 104, 134, 164, 184, 202, 212, 232, 244 and 248; and the largest values are 
C(210) = 19, C(-ff) = 18 for H = 30,60, 102 and 156, ((150) = 16 and ((198) = 13. 
Table 4 shows the growth of Z{H). 



36 J. M. CHICK 

Table 4. Cumulative total of CsN's with gcd{p',q') < H 



H 


20 


40 


60 


80 


100 


120 


140 


160 


180 


200 


220 


240 


260 


268 


Z{H) 


44 


98 


162 


204 


263 


334 


390 


450 


491 


531 


578 


646 


687 


705 



Table 5 shows ^{n : n is a n = c (mod m) and n < 10^} for various m, c 
and N. 

Table 5. Cumulative totals of C^N's up to 10^ satisfying n = c (mod m) 



N 



9 11 



13 



15 



17 



19 



21 



23 



m 


c 




5 


1 


35 


133 


457 


1405 


4611 


14716 


49030 


169157 


595168 




2 


1 


6 


40 


133 


455 


1522 


5151 


17479 


61711 




3 


5 


11 


41 


138 


434 


1421 


4726 


16108 


56953 




4 


3 


19 


49 


179 


580 


1877 


6399 


22016 


77050 


7 


1 


22 


102 


339 


1078 


3472 


11029 


36668 


125774 


443797 




2 


4 


9 


36 


136 


499 


1660 


5590 


19280 


68227 




3 


4 


18 


54 


171 


501 


1636 


5645 


19551 


68150 




4 


3 


12 


55 


162 


544 


1766 


6057 


20990 


73529 




5 


4 


10 


46 


133 


494 


1666 


5379 


18752 


65616 




6 


4 


15 


54 


172 


567 


1776 


5964 


20410 


71560 


11 


1 


13 


48 


183 


591 


2063 


6678 


22417 


77368 


272654 




2 


3 


18 


54 


161 


471 


1499 


4965 


17230 


60546 




3 


2 


13 


43 


134 


432 


1367 


4729 


16599 


58510 




4 


4 


13 


40 


142 


421 


1362 


4598 


15971 


55563 




5 


4 


13 


47 


151 


435 


1411 


4756 


16204 


57647 




6 


6 


15 


49 


151 


478 


1515 


4944 


16670 


58038 




7 


2 


10 


34 


127 


443 


1378 


4711 


16318 


57489 




8 


5 


14 


46 


130 


455 


1511 


4869 


16543 


57988 




9 


4 


12 


48 


115 


420 


1431 


4748 


16233 


57508 




10 


3 


15 


45 


155 


464 


1386 


4571 


15626 


54941 


12 


1 


38 


145 


516 


1632 


5353 


17221 


57694 


199002 


700227 




5 


4 


11 


23 


72 


242 


748 


2456 


8444 


29527 




7 


4 


15 


50 


153 


478 


1517 


4994 


16766 


59215 




11 














9 


52 


164 


550 


1915 



In [12] Pinch describes and searches for certain special types of CiV discussed 
by other authors, including strong Fibonacci pscudoprimes (of which he finds just 
one up to 10^*, with d = 8). These special CN^s all have the property that Pi + 1 
divides n± 1 for I < i < d, and for c? = 3 I have proved that no such numbers exist: 
see [H]. 
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